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ABSTRACT 

We present photometry of the large scale environments of a sample of twelve broad line AGN with 
0.06 < z < 0.37 from deep images in the SDSS u, g, r, and i filters taken with the 90Prime prime 
focus camera on the Steward Observatory Bok Telescope. We measure galaxy clustering around 
these AGN using two standard techniques: correlation amplitude (B g9 ) and the two point correlation 
function. We find average correlation amplitudes for the 10 radio quiet objects in the sample equal 
to (9±18, 144±114, -39±56, 295±260) Mpc 1,7 ' in (it, g, r, i), all consistent with the expectation 
from galaxy clustering. Using a ratio of the galaxy-quasar cross-correlation function to the galaxy 
autocorrelation function, we calculate the relative bias of galaxies and AGN, b gq . The bias in the 
u band, b gq = 3.08 ± 0.51 is larger compared to that calculated in the other bands, but it does not 
correlate with AGN luminosity, black hole mass, or AGN activity via the luminosity of the [O III] 
emission line. Thus ongoing nuclear accretion activity is not reflected in the large scale environments 
from ~10 /i -1 kpc to ^0.5 h~ l Mpc and may indicate a non-merger mode of AGN activity and/or 
a significant delay between galaxy mergers and nuclear activity in this sample of mostly radio quiet 
quasars. 

Subject headings: galaxies: active;large-scale structure of universe 


1. INTRODUCTION 

In local galaxies, the masses of the nuclear 
black holes are correlated with the stellar veloc¬ 
ity dispersion of the host galaxy bulges, the Mrh- 
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Theoretical modeling establishes a framework for the 
coevolution of galaxies and supermassive black holes 
such that after the initial formation of central black 
holes in galaxies, nuclear activity a nd star f ormatio n 
are triggered by a merger event jBarnes_^_JIerMuisd 
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ISomerville et al.ll2008f) . These hierarchical models make 
specific predictions for the large scale environments of 
luminous quasars. 

Observational studies have investigated the wide pa¬ 
rameter space of this problem by examining: the quasar- 
quasar and galaxy-quasar clustering in terms of correla¬ 
tion length and amplitude; nearest neighbor and other 
statistics for estimating environment density; properties 
of galaxies in AGN environments; and the properties of 
AGN host galaxies, covering also a range in redshift and 
intrinsic AGN properties. Many of these studies have 
used Sloan Digital Sky Survey (SDSS) data, thus con¬ 
centrating on the z < 0.3 universe. 

From the observational results, two points of consensus 
have emerged that generally agree with the models: (1) 
AGN activity is associated with galaxy mergers and 
enha nc ed star formation i n the ir hosts dSanders et al.1 
1988t Kauff mann et all 12004 iGuyon et all 20061 


Canalizo et al.1 20071 : Bennert et al.1 2008l~ Coil et al l 

20091 : iVeilleux et al.l l2009t iTeng fe Veilleuxl 1201011 . ah 


though there is increasing evidence, e.g. due to lack 
of correlation of galaxy str uctural parameter s with 
local environment density, (|Kauffmann et, all l2004fl 
and an absence of increased rate of di sturbed host 
morphologies in AGN hos ts up to z ~ 1.3. (|Grogin et al.l 
12004 iGabor et al.1 I2009H . that the onset of nuclear 
activity is delayed with re s pect t o the initial galaxy 
interaction (ITJ et all 120061 12001 lEllison et al.lMTl 
and that there is a significant n on-n r erger mode 
to AGN activity (1 Hopkins fe Hernauistl 1200611 : (2) 

on megap arsec scales, bright AGN cluster like L* 
galaxies (lAdelberger fc Steidell [20051: [IT et al.l 120061 


Serber et al.1 2004 Coldwell fc Lambasl 120061: iCoil et al l 
2007HGilli et al.ll2009l l. and over time, they occupy dark 
matter halos of constant mass (~ 2 —3 x 10 12 h~ x Mg) so 
they were more biased tracers of mass at higher redshift 
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than thev are in the current enoch (Porciani et al. 

2004 

Wake et al.l 12004 

Croom et al. 2005: Mvers et al.1 

2006 

2007al Shen et al. 

2007; Ross et al. 2003) ■ 


pendence on quasar luminosity (Croor 

a et al. 2005 

Porciani & Norberd !2006t Mvers et al. 

2006, 2007a; 

Shen et al.1 

2009 

except at the lowest redshifts 


2004 [Constantin fe Vogelevll2006D . Galaxy 


clustering around AGN does correlate with nuclear ac¬ 
tivity especially for luminous, strongly accreting AGN, 
which tend to lie in massive ho sts (M* > 3 x lO^Mq ) 
with higher than average SFR (jKauffmann et all 1200 3). 
and are more likely to be found i n lower density environ¬ 
ments (i Kauffmann et al .1 12004 iConstantin et al.l 120081 : 
ISilverman et al.1 120091 : ISabater et al.l 120121 I2013T) . No 
quasar luminosity dependence on galaxy-quasar clu ster¬ 
ing is seen in the DEEP2 sample with z ~ 1 ijGoil et alJ 
[2007 1. though quasars are found to cluster like blue, star¬ 
forming galaxies. 

AGN fuelin g via a merger-independen t quiescent or 
Seyfert mode dHopkins fe Hernauistl 120061) may be most 
important for weaker AGN, as they show no prefer- 
ence for low density r egions in low redshift samples 
(|Kauffmann et afl[2004f ) and X-ray selected AGN at z ~ 
1 with lower mass hosts (M* < 10 11 M^1 which are foun d 
in a broad range of environments (ISilverman et al.ll2009f ). 
However, there is evidence that the transit ion from 
Seyferts to LINERS proposed by iKewlev et al.l (1200611 is 
affected by the dens ity of the environment, in SDSS sam¬ 
ples at low redshift (|Constantin et ahll2008 ]_ and inthc X- 
ray se lected z ~ 1 AGN s ample of MrmterozDortaeLalJ 
(12009( 1. These authors. and IGonstantin fc Vogelevl ( 2006( 1 
also Hnd that that the LINERS cluster more strongly 
than Seyfert galaxies. 

Unlike galaxies, quasars do show excess clu stering on 
small scales (10-100 ft." 1 kpc) up to z ~ 3 dHennawi et al.l 
120061 iMvers et al.l l2007bL IShen et all 1 201 011 ~ And while 
studies of the ~1 ft . -1 Mpc environments around low 
redshift AGN show no distinct correlation between nu¬ 
clear activity a nd th e presence of galaxy neighbors 
on these scales (fSerber et al.l 120061 : iColdwell fc Lambasl 
120061 iLi et al.l I2008H . there is evidence that the pres¬ 
ence of a galaxy comp anion on smaller s cales does in¬ 
fluen ce AGN activity (| Alonso et al.l 120071 : ISabater et al.1 
I2013D . On larger scales, z > 1 radio loud quasars 
are fo und in richer en v ironments than rad i o quiet ob¬ 
jects (iHall et all 119981 iTeplitz et all 119991 iWold et, all 
2003) while for low redshift this distinction is less 
pronounced and AGN generally are not p referent ially 
found at the centers of r i ch clusters dFisher et all 
19961 iCroom fe Shanks! 119991 IMcLure fe Dunlopl 2001b : 


Wold et all (20011 Goldwell et al.l 120028 : iBarr et all 12003 : 

Goldwell fc Lambasl 120061 1 . For cold. or QSO, mode ac¬ 


cretion, this can be understood as inhibition due to 
the stripping and harassment that occurs in these en¬ 
vironments. However, some studies have shown that 


quasars do trace the large s c ale st r ucture of galaxy 
clusters (|Sochting et al.l 120021 12004 iGilli et aD 120031 


ISilverman et al.ll2008H on scales up to 10 h -1 Mpc. 

In this paper, we investigate the large scale environ¬ 
ments of a sample of 12 relatively bright low redshift 
(0.06 < z < 0.37) broad line AGN using deep images in 
the SDSS u, g , r, and i filters. The fields presented here 


overlap almost entirely with the SDSS, but our photome¬ 
try goes significantly deeper, in some fields and filters by 
up to 2.5-3 magnitudes. Our wide field multiband study 
on scales ~1 degree, is distinct from previous single filter 
studies with Hubble Space Telescope (HST) using the ^3 
areminute field of t he Wide Field Planetary Camera 2 
(I Fisher et al.1 119961 : iFinn et all 120011 : IMcLure &: Dunlopl 
I2001bl) . Likewise, the multiband nature of our study 
extends previous wider field studies usi ng t he Bj and 
R data from UK Schmidt plates dBrown et al.l 120011 : 
IColdwell et al.l[20021: ISochting et al.ll2002h . This study is 
a blend of a statistical approach and a detailed study of 
each field. We calculate the correlation lengths and am¬ 
plitudes of the galaxy-quasar clustering and compare re¬ 
sults in different bands. We also use our data and results 
from the literature to discuss each quasar’s environment 
in detail. Our primary aim is to probe galaxies as far 
down the luminosity function as allowed by the depth of 
the photometry to as wide an area around each quasar 
as possible in order to test if galaxy clustering around 
quasars is comparable to that of galaxy-galaxy clustering 
to these depths at these scales. Throughout this paper 
we assume a cosmology with h = 0.71 (Hq = 100 h km 
s -1 Mpc -1 ), = 0.27, and Qa = 0.73. 

2. DATA 

2.1. Observations and Reductions 

We made our observations with the 90Prime prime fo¬ 
cus imager (IWilliams et al.1 (2004 1 on the 2.3-meter Bok 
Telescope at Kitt Peak National Observatory in January 
- June, 2008. 90Prime is a wide field camera consist¬ 
ing of four CCDs arranged in a square mosaic located at 
the prime focus of the telescope with a field of view of 
~1 degree. We used the SDSS it, g, r, and i filters for 
the photometry, and with the exception of HS0624+690, 
all our fields lie inside the SDSS footprint. However, our 
photometry is significantly deeper than that of the SDSS 
survey, as discussed further in the next section. We chose 
the quasars in our sample to have galaxy environments 
accessible to deep imaging with our 2-meter class tele¬ 
scope and with ultraviolet spectra in the archives of the 
Hubble Space Telescope or Far Ultraviolet Spectroscopic 
Explorer for future studies of line of sight and associ¬ 
ated absorption in these systems. The lower bound on 
the redshifts of the sample quasars was determined by 
the size of the chips in the CCD mosaic and the ~Mpc 
scales targeted in this study. We compiled a sample of 
~20 target fields for the first semester of the 2008 observ¬ 
ing season and completed the photometry in at least two 
bands for 12 of those fields. Total exposure times were 
estimated for each field and filter in order to achieve sig¬ 
nificant detections of galaxies with M* + 2 at redshift of 
the quasar in each filter. The final sample presented here 
includes one Seyfert, 9 radio quiet quasars and two radio 
loud quasars and so somewhat overrepresents the radio 
loud population, although we separate these fields out in 
the discussion of our results. Follow up observations to 
obtain data on other fields of quasars with available UV 
spectra and with appropriate redshifts are possible. 

The quasars themselves were centered on one of the 
chips of the mosaic, usually chip 1, rather than in the 
center of the field due the 8.3' interchip gaps and several 
traps and other quality issues with the other chips in the 
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mosaic at the time of the observations. The raw data 
frames are multiextension FITS files consisting of eight 
2048x4096 images, one for each of the two amplifiers on 
each chip. Amplifiers 1 and 2 lie on chip 1, 3 and 4 on 
chip 2, 5 and 6 on chip 3, and amps 7 and 8 are on chip 
4. We give a summary of the observations in Table [l] 
We reduced the data with standard techniques for bias, 
zero, and flat field corrections, including fringe correc¬ 
tions in the i band, using the standard reduction tasks 
in the IRAF mosaic package. 




TABLE 1 

90Prime Observations 


QSO 

Type 1 

RA 

Dec 

Redshift 

r 

u exp.(s) 

9 exp.(s) 

r exp.(s) 

i exp.(s) 

MRK586 

Syl.2 

02:07:49.8 

+02:42:55 

0.156 

15.42 

2800 

4400 

3600 


HS0624+690 

QSO 

06:30:02.5 

+69:05:04 

0.370 

14.73 2 

9200 

6400 

4400 

1200 

PG0844+349 

Syl.O 

08:47:42.4 

+34:45:04 

0.064 

14.48 

4400 


7600 


PG0923+201 

Syl.O 

09:25:54.7 

+19:54:05 

0.192 

15.55 

5600 

5200 

2000 

2400 

PG0953+414 

Syl.O 

09:56:52.3 

+41:15:22 

0.234 

14.93 

8000 

2400 

3600 

4200 

PG1116+215 

Syl.O 

11:19:08.7 

+21:19:18 

0.177 

14.46 

5370 

400 

400 

2400 

PG1307+085 

Syl.2 

13:09:47.0 

+08:19:48 

0.155 

15.59 

4400 

2000 

3200 

2400 

PG1404+226 

NLS1 

14:06:21.9 

+22:23:46 

0.098 

16.02 

5200 

4800 

6800 

2400 

PG1444+407 

Syl.O 

14:46:45.9 

+40:35:06 

0.267 

15.81 

4400 

4000 

3200 

2400 

PG1545+210 

Syl.2 

15:47:43.5 

+20:52:17 

0.264 

15.56 

6000 

4800 

2800 

2400 

PG1612+261 

Syl.5 

16:14:13.2 

+26:04:16 

0.131 

15.87 

4000 

4800 

3800 

2400 

Q2141+175 

Syl.O 

21:43:35.5 

+17:43:49 

0.211 

16.04 


4500 

3000 

1500 


1 IVeron-Cettv &; Veronl (|2010f) 
^ Measured from 90Prime data 
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Fig. 1.— Example 90Prime data frame. The PG0923+201 field 
in u,g (upper left, right) and r,i (lower left,right). 


2.2. Photometric Calibration and Catalog 
Construction 

2.2.1. Calibration and Masking 

We flux calibrated each image to the SDSS images of 
the same fields by selecting stars with magnitudes fainter 
than ~ 15 to avoid saturation, and brighter than ~20 
where the asinh magnitud e use d by SDS S diverges from 
the traditional magnitude (ILupton et al.ll 19991 ). We used 
this sample to estimate the zero point magnitudes and 
the first order extinction coefficients for each 90Prime 
field and then applied these to the data. We combined 
the flux-calibrated images of each field to create a full 
exposu re-t im e co- added image for each filter using SWarp 
(|Bertin et alJ 1200 2). We cross-checked star and galaxy 
magnitudes once again with SDSS galaxies with r < 21 to 
verify the final quality of the calibration in the co-added 
images. Calibrated u,g,r, and i frames for PG0923+201 
are shown in Fig. [1] 

HS0624+6907 does not overlap with the SDSS foot¬ 
print, so to calibrate the frames of this field, we use the 
frames taken in each filter on photometric nights as ref¬ 
erences for the frames taken on non-photometric nights. 
The zero point magnitude and first order extinction co¬ 
efficient is calculated from the SDSS fields taken on the 
same photometric nights, and so this photometry is also 
tied to the SDSS system as the rest of the sample. We 
calculated photometric, Kron, a nd Petrosian magnitudes 
using the SExtractor package (IBertin fe Arnoutslll996D 
on the full exposure-time flux-calibrated co-added im¬ 
ages of each filter. We found that optimal source extrac¬ 
tion was achieved by considering each amplifier section of 
each image on a case-by-base basis for amps 1-8 and the 
extraction parameters used thus vary slightly over the 
full mosaic. We give a brief summary of those parameter 
choices here. Amps 1, 2, and 3, are typically low noise 
and free of bad columns. So we used 1.5cr for detection 
and analyze thresholds, DETECT_THRESH and ANA- 
LYZE_THRESH. Since the fields are crowded, we use 
a balanced combination of background and deblending 
parameters, BACK.SIZE, DEBLEND_NTHRESH, and 
DEBLEND_MINCONT, to detect as many individual 
objects and as close as possible to stars or foreground 
galaxies but to avoid over-resolving. Amp 4 has noise 
patterns which are enhanced by resampling. To over¬ 



Fig. 2.— Example of a frame mask, for the PG0923+201 in u 
filter, shown in Fig. [T] including inter-chip gaps, trap in amps 5 
and 6, bad pixel/columns, and bright stars. 


come this problem we used a moderately higher detec¬ 
tion threshold, ~3, and larger area to estimate the back¬ 
ground with deblending parameters set to identify the 
strips as artificial bright sources. In some extreme cases 
we mask the noise pattern manually. Amps 5 and 6 have 
a higher noise level than amps 1 and 2 but typically free 
of fixed patterns, so we use the same technique as amp 
1 but with a higher detection threshold, ~1.7. Amps 7 
and 8 behave generally like amps 5 and 6. We dealt with 
a few cases of noise patterns using the same technique as 
amp 4, with detection a threshold of ~2.3. 

To remove bad columns and bright stars and their as¬ 
sociated diffraction patterns, we created image masks 
which consist of: (1) a manually generated mask for the 
interchip gaps; (2) a mask for a consistent trap present 
in chip 3, overlapping amps 5 and 6; (3) a bright star 
mask created manually by identifying the bright stars 
in a given frame from the SExtractor output and grow¬ 
ing the masked regions using the measured FWHM; (4) a 
bad pixel mask created using the IRAF task obj masks to 
mask pixels with pixel values 15cr below the background; 
and finally, (5) a bleed trail mask created by using the 
segmentation map created by SExtractor for a detection 
threshold ^80-100 and very low debelending parameter 
and large background estimates, thus including all the 
bleed patterns. We created masks for each image of each 
field in every filter individually. The downside of this 
technique is that step (5) can include some bright fore¬ 
ground galaxies into the mask, so we recover these from 
careful inspection of the images in all four filters. As 
an example, our mask frame for PG0923+201 in u filter 
(Fig. [T]) is shown in Fig. [2] Due to the large dithering 
of individual frames before co-adding, some regions have 
lower co-added exposure time than the full stacked im¬ 
age. To decrease the biases on galaxy counting when 
using a certain detection limit on magnitude, we cut the 
corners and more from inner-gaps to include only the re¬ 
gions with uniform co-added exposure time in each final 
image. We also mask objects with centers within 8 pix¬ 
els of a masked region or an edge from the final catalogs. 
Finally we filter our catalog sources based on the internal 
SExtractor flags. Here, we present the catalog of objects 
detected at 5cr or greater in r and at 3<r or greater in 
it, g and i in Table [2] as a machine readable table. The 
magnitudes fisted are the SExtractor MAG_BEST val¬ 
ues, where the Kron_fact and min_radius parameters gov- 
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erning MAG_AUTO were set to 2.5 and 1.5, respectively. 
Magnitudes of 99 in this table denote non-detections in 
the u, g, or i filter, and values of -99 label objects not 
observed in one of those filters, primarily due to masking. 
In the former case, the error bar gives the lcr detection 
limit for later use with photometric redslrift codes. For 
the analyses discussed in this paper, we use only galaxies 
detected at 5cr significance or greater. 



TABLE 2 

90Prime Catalog of Galaxies in Quasar Fields 


QSO field 

RA 

Dec 

u 1 

&U 

9 


r 

CJ y 

i 


A 2 


B 2 

CB 

radius 3 

9 3 

FWHM 3 

class 3 

HS0624p690 

98.0779 

68.8485 

23.10 

0.14 

99.00 

25.52 

21.20 

0.08 

20.65 

0.11 

1.443 

0.109 

0.840 

0.064 

2.56 

-12.91 

5.05 

0.02 

HS0624p690 

98.1246 

68.8445 

99.00 

27.74 

99.00 

25.52 

22.67 

0.18 

99.00 

25.54 

0.677 

0.143 

0.519 

0.110 

7.79 

-75.60 

2.93 

0.00 

HS0624p690 

97.2567 

68.8409 

99.00 

27.74 

99.00 

25.52 

21.93 

0.11 

21.28 

0.16 

0.852 

0.121 

0.531 

0.075 

5.67 

-85.21 

2.72 

0.61 

HS0624p690 

98.0474 

68.8475 

99.00 

27.74 

99.00 

25.52 

20.80 

0.06 

19.90 

0.08 

1.213 

0.080 

0.863 

0.057 

3.97 

21.11 

2.96 

0.96 

HS0624p690 

97.5498 

68.8449 

99.00 

27.74 

99.00 

25.52 

21.03 

0.07 

99.00 

25.54 

0.964 

0.069 

0.920 

0.065 

5.46 

-35.21 

2.37 

0.26 

HS0624p690 

97.6265 

68.8442 

99.00 

27.74 

99.00 

25.52 

22.11 

0.13 

22.12 

0.23 

0.799 

0.116 

0.603 

0.088 

5.97 

61.78 

2.45 

0.57 

HS0624p690 

97.4941 

68.8537 

99.00 

27.74 

99.00 

25.52 

21.78 

0.11 

21.71 

0.19 

0.981 

0.118 

0.682 

0.082 

5.02 

-82.06 

3.82 

0.12 

HS0624p690 

98.1444 

68.8458 

99.00 

27.74 

99.00 

25.52 

21.85 

0.11 

20.80 

0.12 

0.799 

0.107 

0.562 

0.075 

5.90 

73.74 

2.28 

0.73 

HS0624p690 

97.5849 

68.8437 

99.00 

27.74 

99.00 

25.52 

22.75 

0.17 

22.25 

0.24 

0.535 

0.129 

0.432 

0.104 

6.39 

59.89 

2.28 

0.50 

HS0624p690 

97.7454 

68.8461 

99.00 

27.74 

99.00 

25.52 

21.79 

0.11 

20.74 

0.12 

0.858 

0.109 

0.740 

0.094 

5.77 

29.80 

2.87 

0.10 

HS0624p690 

97.2771 

68.8423 

99.00 

27.74 

99.00 

25.52 

22.81 

0.18 

99.00 

25.54 

0.466 

0.144 

0.402 

0.125 

8.71 

71.16 

2.58 

0.34 

HS0624p690 

97.8299 

68.8482 

99.00 

27.74 

99.00 

25.52 

21.17 

0.08 

20.56 

0.11 

1.254 

0.098 

0.846 

0.066 

3.67 

46.78 

3.74 

0.21 

HS0624p690 

97.5091 

68.8468 

99.00 

27.74 

21.62 

0.10 

21.40 

0.09 

21.00 

0.14 

0.942 

0.086 

0.733 

0.067 

4.76 

-29.61 

2.48 

0.86 

HS0624p690 

97.5560 

68.8461 

99.00 

27.74 

99.00 

25.52 

22.41 

0.15 

99.00 

25.54 

0.642 

0.111 

0.559 

0.097 

6.31 

18.58 

2.16 

0.64 

HS0624p690 

98.0133 

68.8473 

99.00 

27.74 

99.00 

25.52 

22.77 

0.18 

21.12 

0.14 

0.753 

0.166 

0.462 

0.102 

5.87 

-88.26 

3.53 

0.23 

HS0624p690 

97.3695 

68.8452 

99.00 

27.74 

99.00 

25.52 

22.90 

0.18 

22.11 

0.23 

0.704 

0.156 

0.423 

0.094 

5.34 

75.93 

3.26 

0.52 

HS0624p690 

97.3791 

68.8490 

99.00 

27.74 

99.00 

25.52 

20.78 

0.06 

99.00 

25.54 

1.245 

0.079 

0.874 

0.055 

5.03 

28.43 

2.92 

0.20 

HS0624p690 

97.9586 

68.8496 

99.00 

27.74 

99.00 

25.52 

22.51 

0.16 

21.22 

0.15 

0.571 

0.125 

0.490 

0.107 

6.94 

75.39 

2.82 

0.49 

HS0624p690 

97.5345 

68.8568 

19.66 

0.03 

18.32 

0.02 

17.46 

0.01 

17.10 

0.02 

2.327 

0.029 

2.117 

0.027 

3.37 

-74.43 

3.91 

0.03 

HS0624p690 

97.5649 

68.8494 

99.00 

27.74 

99.00 

25.52 

22.12 

0.13 

21.59 

0.18 

0.542 

0.155 

0.378 

0.108 

8.39 

-13.78 

2.36 

0.36 

HS0624p690 

97.0202 

68.8484 

99.00 

27.74 

99.00 

25.52 

20.98 

0.07 

20.17 

0.09 

1.610 

0.121 

0.900 

0.068 

4.66 

54.87 

4.81 

0.39 

HS0624p690 

97.8307 

68.8524 

99.00 

27.74 

99.00 

25.52 

21.14 

0.08 

20.04 

0.09 

0.858 

0.065 

0.760 

0.057 

4.07 

79.73 

2.00 

0.96 

HS0624p690 

97.3351 

68.8497 

22.32 

0.09 

99.00 

25.52 

21.71 

0.10 

21.37 

0.16 

0.953 

0.101 

0.758 

0.080 

3.94 

65.93 

3.44 

0.08 

HS0624p690 

98.1068 

68.8543 

99.00 

27.74 

20.92 

0.07 

19.92 

0.04 

18.59 

0.04 

1.174 

0.047 

0.915 

0.037 

3.56 

80.00 

2.86 

0.94 

HS0624p690 

97.1526 

68.8478 

99.00 

27.74 

99.00 

25.52 

22.41 

0.15 

99.00 

25.54 

0.676 

0.147 

0.501 

0.109 

7.64 

-77.44 

3.48 

0.47 

HS0624p690 

97.9007 

68.8552 

99.00 

27.74 

99.00 

25.52 

20.11 

0.05 

19.11 

0.06 

1.381 

0.062 

0.896 

0.040 

3.47 

9.03 

3.23 

0.94 

HS0624p690 

97.7168 

68.8546 

99.00 

27.74 

99.00 

25.52 

20.83 

0.07 

20.13 

0.09 

1.116 

0.072 

0.897 

0.058 

4.24 

-75.39 

2.69 

0.53 

HS0624p690 

97.8514 

68.8552 

99.00 

27.74 

99.00 

25.52 

21.59 

0.10 

21.05 

0.14 

0.912 

0.125 

0.780 

0.107 

6.80 

0.23 

4.54 

0.00 

HS0624p690 

97.9116 

68.8545 

99.00 

27.74 

99.00 

25.52 

22.63 

0.16 

21.52 

0.17 

0.705 

0.127 

0.522 

0.094 

4.74 

74.68 

3.43 

0.64 

HS0624p690 

98.0572 

68.8560 

22.87 

0.13 

99.00 

25.52 

21.29 

0.09 

20.81 

0.12 

1.592 

0.142 

0.783 

0.070 

4.98 

87.82 

5.97 

0.00 

HS0624p690 

97.6051 

68.8607 

99.00 

27.74 

99.00 

25.52 

22.02 

0.12 

22.28 

0.24 

0.863 

0.114 

0.664 

0.088 

5.13 

85.51 

2.90 

0.12 


1 MAG_BEST from SExtractor 
^ semimajor and semiminor axes in arcsec 

3 SExtractor parameters KRON_RADIUS, THETA_J2000 (in degrees), FWHM, and CLASS-STAR as discussed in the text and in lBertin et all J2002D 
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TABLE 3 

Flux Calibration Errors 


TABLE 4 

Individual Flux Calibration Errors 


filter 

Systematic err. 

95% conf. 

Total 

95% conf. 

u 

-0.039 

0.011 

0.305 

0.009 

g 

-0.083 

0.006 

0.152 

0.005 

r 

-0.028 

0.008 

0.122 

0.007 

i 

-0.022 

0.013 

0.066 

0.006 



Model Mag 90Prime 

Fig. 3.— Example of error estimates from photometric calibra¬ 
tion process in g filter from all fields. 

We estimate systematic error and scatter due photo¬ 
metric calibration for all four filters using a final compar¬ 
ison with the SDSS photometry. For this comparison, we 
use all objects detected in our catalogs at >5er and with a 
SExtractor classification parameter (CLASS.STAR, de¬ 
scribed in the next section) likely to be a galaxy, <0.3. 
We also require the relative error in the SDSS magni¬ 
tude to be smaller than 1%. The values of the systematic 
shifts are small, typically less than 0.04. They are listed 
in Table [3] and plotted for the g band in Fig. [3] Also 
in Table [3] we list the 95% confidence on the systematic 
offset, along with the total photometric statistical un¬ 
certainties in each band and their 95% confidence inter¬ 
val. The total statistical uncertainties are obtained from 
adding the systematic offsets from SDSS in quadrature 
with the standard deviation of the distribution of SEx¬ 
tractor magnitude errors for each filter. The mean and 
standard deviation of those distributions are listed sepa¬ 
rately for each field in Table 0] 

2.2.2. Star-Galaxy Separation and Filtering 

SExtractor calculates a CLASS_STAR flag based on 
core geometry of the light distribution of the detected 
objects, where values near 1 indicate a centrally concen¬ 
trated star and those near 0 indicate a galaxy. This pa¬ 
rameter is sensitive to the value of the input parameter 
SEEING_FWHM. We measure this input parameter from 
objects clearly classified as stars (CLASS-STAR > 0.9) 
in initial photometric catalogs, which were created from a 
first pass with SExtractor on the masked images. To im¬ 
prove the accuracy of the CLASS-STAR star-galaxy sep¬ 
aration output parameter, we performed this measure¬ 
ment iteratively until we reached convergence between 
the input value of the SEEING-FWHM and the mea¬ 
sured FWHM of the stellar catalog sources. We achieved 
convergence of <0.1" in 1-2 iterations. In order to deter¬ 


field 

filter 

Systematic err. 

95% conf. 

Total 

95% conf. 

MRK586 

u 

-0.061 

0.016 

0.328 

0.012 

MR.K586 

g 

-0.026 

0.012 

0.253 

0.010 

MRK586 

r 

0.007 

0.016 

0.139 

0.014 

PG0844+349 

u 

0.009 

0.013 

0.293 

0.010 

PG0844+349 

r 

-0.031 

0.008 

0.131 

0.007 

PG0923+201 

u 

0.003 

0.013 

0.356 

0.011 

PG0923+201 

g 

-0.062 

0.012 

0.169 

0.011 

PG0923+201 

r 

0.022 

0.012 

0.144 

0.010 

PG0923+201 

i 

0.119 

0.020 

0.199 

0.016 

PG0953+414 

u 

0.005 

0.011 

0.286 

0.009 

PG0953+414 

g 

-0.030 

0.008 

0.134 

0.006 

PG0953+414 

r 

-0.008 

0.010 

0.114 

0.009 

PG0953+414 

i 

0.062 

0.013 

0.131 

0.010 

PG1116+215 

u 

-0.019 

0.013 

0.317 

0.010 

PG1116+215 

g 

-0.029 

0.025 

0.118 

0.015 

PG1116+215 

r 

-0.019 

0.041 

0.070 

0.019 

PG1116+215 

i 

0.095 

0.018 

0.171 

0.014 

PG1307+085 

u 

0.034 

0.010 

0.292 

0.008 

PG1307+085 

g 

-0.046 

0.010 

0.148 

0.008 

PG1307+085 

r 

-0.014 

0.006 

0.071 

0.005 

PG1307+085 

i 

0.067 

0.018 

0.145 

0.012 

PG1404+226 

u 

-0.117 

0.016 

0.376 

0.013 

PG1404+226 

g 

-0.076 

0.010 

0.172 

0.008 

PG1404+226 

r 

-0.032 

0.008 

0.134 

0.007 

PG1404+226 

i 

0.064 

0.022 

0.171 

0.018 

PG1444+407 

u 

-0.083 

0.014 

0.284 

0.012 

PG1444+407 

g 

-0.089 

0.008 

0.148 

0.007 

PG1444+407 

r 

-0.036 

0.018 

0.111 

0.016 

PG1444+407 

i 

0.006 

0.013 

0.145 

0.011 

PG1545+210 

u 

-0.005 

0.007 

0.331 

0.006 

PG1545+210 

g 

-0.091 

0.014 

0.165 

0.013 

PG1545+210 

r 

-0.017 

0.011 

0.093 

0.008 

PG1545+210 

i 

0.036 

0.017 

0.161 

0.014 

PG1612+261 

u 

-0.088 

0.014 

0.305 

0.011 

PG1612+261 

g 

-0.110 

0.010 

0.128 

0.008 

PG1612+261 

r 

-0.045 

0.011 

0.108 

0.008 

PG1612+261 

i 

0.033 

0.012 

0.185 

0.010 

Q2141+175 

g 

-0.093 

0.010 

0.155 

0.007 

Q2141+175 

r 

-0.058 

0.012 

0.100 

0.011 

Q2141+175 

i 

0.010 

0.021 

0.185 

0.014 


mine how to use this CLASS-STAR parameter to best 
distinguish stars from galaxies in our final catalogs, we 
performed Monte-Carlo simulations of our data set using 
synthetic input stars that we then extract from the data 
frames using the same SExtractor parameters as we used 
on the data. Using the IRAF task mkobjects, we gen¬ 
erated synthetic stars with magnitudes between 15 and 
27 and with PSFs that matched those of stellar objects 
in the post-iteration catalogs. We randomly distributed 
these synthetic stars on chips 1 and 4 (amps 1,2,7, and 
8) of each image, to bracket the best and worst quality 
CCDs in the mosaic. After extracting these simulated 
stars, we see the typical behavior of the CLASS-STAR 
parameter, i. e. that it is a reliable star-galaxy separa¬ 
tor at bright magnitudes but at a particular threshold 
magnitude, CLASS-STAR is spread evenly between 1.0 
and some saturation value, 0.35 for our 90Prime frames. 
Thus, at magnitudes dimmer than the threshold, we can¬ 
not rely on this parameter alone to separate galaxies from 
stars. This threshold magnitude varies across our data 
frames, but for each frame, we fit a function to the lower 
envelope of CLASS-STAR versus magnitude, m: 

CLASS-STAR = a/(1.0 + exp((m - b)/c)) + d (1) 
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i 


< 
h- 
C O 

cn 0.6 


0.2 

0 

Fig. 4.— CLASS_STAR versus r for PG1612+261 with the fit 
to Equ. [T] Cyan(red) points are points for which the difference 
between the input and output stellar magnitudes in the Monte- 
Carlo simulations was less than 0.8(0.2). The fit is designed to 
find the best inner envelope for the red data reaching asymptotic 
value of 0.95 at lowest magnitude. 

where m is the magnitude in the relevant filter, and a, 
b, c, and d are fitted parameters. We use the fits to this 
expression in our final star-galaxy separation algorithm. 
Figure |4] shows an example of this fitting. 

Equ. [I] is used to define a cut in the CLASS-STAR 
parameter as function of source magnitude. We find 
that using this function as a strict cut results in many 
faint sources lost from our catalogs, so we relax this 
criterion and rely on a color based criterion as well. 
To construct the color criterion, we use the unambigu¬ 
ous stellar sources in the matched catalogs, those with 
CLASS-STAR > 0.95 in all filters, along with the sources 
in the Gunn-Stryker Atlas (iGunn fe Strvkedfl983D to de¬ 
fine the stellar locus in color-color space. We then ex¬ 
clude from our catalogs any source that lies within a 
specified distance from that locus in both g — r versus 
u — g and r — i versus g — r. This distance is allowed to 
vary in order to give the best match in number counts 
to our adopted luminosity functions as described below, 
but it is typically 0.1-0.3. In order to use the color in¬ 
formation for sources with only one availabl e co l or, we 
use a technique similar to that described by [Coil et all 
1 2004 ). The binned distributions of the stars in our cat¬ 
alogs in each color, normalized to have a maximum value 
of unity, are used to assign a probability of being a star, 
p s . Objects having a p s above a threshold value are ex¬ 
cluded from the final galaxy catalogs. This threshold is 
also allowed to vary slightly, typically between 0.7 and 
0.8, to get good agreement in the number counts in each 
filter. The resulting number counts from these final cat¬ 
alogs, shown below, give us confidence in this star-galaxy 
separation technique. In Figure [5j we show the positions 
of the 5 a galaxies in u, g , r, and i in the portion of the 
field of PG1444+407 covered by chip 1, along with the 
position of the quasar itself. 

2.2.3. Limiting Magnitudes 

We characterize the depth of our data frames in two 
ways: (1) by finding the magnitude at which the lower 
envelope of the distribution of magnitude errors for stars 
in the catalogs equals 0.217, or 5a and (2) by using the 
magnitude distribution of all catalog objects detected at 


8 m <0.80 
8 m <0.20 
f =a/(1.0+exp((m-b)/c))+cF 




6 7 8 9 

X (1000 pix) 


Fig. 5. — Blue, green, red, and orange circles show the positions 
of galaxies detected in u, g, r, and i at > 5cr significance with 
magnitude brighter than the limiting magnitude listed in column 
3 of Table[5]for chip 1 of the PG1444+407 field. Black cross marks 
the position of PG1444+407 itself. 



Model mag 

Fig. 6.— Example of method (1) for calculating limiting magni¬ 
tudes at various significance levels. 


> 4.5cr. We show examples of each in Figures [6] and 
□ and tabulate these values for each co-added frame in 
Table [5] Because we used slightly modified SExtractor 
parameters for the different amps of the full mosaic, we 
calculated the limiting magnitudes separately for each 
amp using method (1). We find that the values for all 
other amps agree with that of amp 1 to within < 0.03 
magnitudes on average, so we report only values for that 
amplifier. 

We treat method (1) as our estimator of the limiting 
magnitude at a particular level of significance. Because 
we use stars only for method (1), it is affected by our 
ability to recover faint extended galaxies from the images. 
Method (2) therefore best characterizes the depth of the 
galaxy catalogs, since those objects dominate the number 
counts at the faintest magnitudes, assuming our star- 
galaxy separation technique is robust. We discuss this 
further below, in Section [3] 

Results for each field/filter with each method are listed 
in Table [5] For reference, the SDSS limiting magnitudes 
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TABLE 5 

Limiting Magnitudes 


QSO 

filter 

5(7 mum 
method (1) 

3(7 miim 
method(l) 

m (complete) 
method (2) 

HS0624+6907 

u 

23.8 

23.9 

22.9 

HS0624+6907 

9 

21.9 

21.9 

20.7 

HS0624+6907 

r 

23.1 

23.4 

22.1 

HS0624+6907 

i 

22.0 

23.1 

21.7 

MRK586 

u 

22.9 

24.0 

22.5 

MRK586 

9 

23.6 

24.5 

22.9 

MRK586 

r 

23.3 

24.4 

23.1 

PG0844+349 

u 

23.9 

24.5 

23.3 

PG0844+349 

r 

24.0 

25.1 

23.3 

PG0923+201 

u 

23.7 

24.8 

23.5 

PG0923+201 

9 

23.6 

24.7 

23.3 

PG0923+201 

r 

22.7 

23.8 

22.5 

PG0923+201 

i 

22.8 

24.0 

22.3 

PG0953+414 

u 

24.1 

25.2 

23.9 

PG0953+414 

9 

22.9 

24.1 

22.9 

PG0953+414 

r 

23.3 

24.4 

23.1 

PG0953+414 

i 

23.6 

24.7 

22.7 

PG1116+215 

u 

23.7 

24.8 

23.5 

PG1116+215 

9 

22.0 

23.0 

21.7 

PG1116+215 

r 

22.2 

23.3 

22.1 

PG1116+215 

i 

22.8 

23.9 

22.5 

PG1307+085 

u 

23.5 

24.6 

23.3 

PG1307+085 

9 

22.7 

23.9 

22.7 

PG1307+085 

r 

23.1 

24.2 

23.1 

PG1307+085 

i 

22.9 

24.0 

22.7 

PG1404+226 

u 

23.6 

24.6 

23.3 

PG1404+226 

9 

23.6 

24.8 

23.5 

PG1404+226 

r 

24.0 

25.1 

23.5 

PG1404+226 

i 

22.8 

23.9 

22.5 

PG1444+407 

u 

23.4 

24.6 

23.3 

PG1444+407 

9 

23.4 

24.6 

23.3 

PG1444+407 

r 

23.2 

24.4 

23.1 

PG1444+407 

i 

22.8 

24.0 

22.7 

PG1545+210 

u 

23.8 

24.9 

23.5 

PG1545+210 

9 

23.5 

24.7 

23.3 

PG1545+210 

r 

23.1 

24.2 

22.9 

PG1545+210 

i 

22.8 

23.9 

22.3 

PG1612+261 

u 

23.3 

24.4 

23.1 

PG1612+261 

9 

23.9 

25.0 

23.7 

PG1612+261 

r 

23.6 

24.7 

23.5 

PG1612+261 

i 

23.1 

24.3 

22.9 

Q2141+175 

9 

23.8 

24.9 

23.7 

Q2141+175 

r 

23.4 

24.5 

23.3 

Q2141+175 

i 

22.6 

23.8 

22.3 


in ( u,g,r,i ) are (22.0, 22.2, 22.2, 21.3), corresponding to 
95% completeness limits. These values may be compared 
with column 4, where we list the 3 a limiting magnitudes 
found via method (1). The stacked g band image of the 
HS0624+6907 field was created from a combination of 
data taken on photometric and non-photometric nights, 
calibrated as described above. The resuling image in this 
particular field and filter had a resulting level of noise 
that necessitated a high detection threshold within SEx- 
tractor to avoid numerous spurious detections around 
bright sources. Thus, the 3- and 5 a limiting magnitudes 
are reported here to be the same, as there are no sources 
extracted at less than 4<r. This also results in a shallow 
field compared to the other photometry reported here. 
Note that the limiting magnitudes in the 90Prime data 
are significantly fainter than SDSS in most other fields 
and filters. 


3. NUMBER COUNTS 

In Figure [8] we show the galaxy number counts in all 
four bins from our 90Prime fields, over a total of 7.35, 
7.45, 8.95, and 7.03 deg 2 in u, g , r, and i, respectively. 



Fig. 7.— Example of method (3) for calculating completeness 
using the magnitude distribution of all catalog objects detected at 
> 4.5 <t. 



Fig. 8. — Number of galaxies in 0.5 magnitude bins per square 
degree in our 90Prime fields in u , g , r, and i. Solid lines are the 
number counts predicted by the luminosity function in each band. 
The vertical dotted and dashed lines show the limiting magnitudes 
for Samples 1 and 2 respectively, listed in Table [6] and discussed in 
Section [4] 


The solid lines in Figure [H] show the predicted num¬ 
ber of galaxies for our adopted luminosity functions in 
each band. Luminosity function parameters adopted are 
shown in Fig. [9] The curves shown here are polyno¬ 
mial fits to the parameter values derived as a function of 
redshift from se vera l sur veys : i n u, w e adopt the param¬ 
eters found by iPrescott et all (|2009f) from the DEEP2 
and SDSS u-band Galaxy Survey' and for the r band, we 
adopt the parameters found bv iCool et all ( 2012 ) from 
the AGES survey. For the g band, we adopt t he B band 
values from the DEEP2 survey (|Faber et all 120071 ). us¬ 
ing the q band lumi nosi ty function parameter measure¬ 
ments of iLovedav et al.1 12012T ) to constrain our fits at 
z = 0. We make this choice because: (a) the wavelength 
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TABLE 6 
Sample Limiting 
Magnitudes 


11 



0 0.5 1 1.5 2 

Redshift 

Fig. 9.— Polynomial fits (solid curves) to lum inosity fun ction pa¬ 
rameter valu es <f>* and M* (squ ares) from Prescott et al. ( 20 03) (u 
band, blue), ([Faber et al .120071) (. B band, gre en).[C ool et al.l (|2012|) 
(r band, red), and For i band (orange) the [Loveday et al. (2012) 
parametrized model for M* is set at the z = 0.6 value for larger 
redshifts, and for 0* it follows the redshift e volution of the r band 
from its z = 0 value. Dashed curves show I Lovedav et al.l (120121) 
parametrized models for the redshift evolution of each parameter 
in i extrapolated to z > 0.5. 


coverage of B overlaps that of the SDSS g filter; (b) 
the agreement between these B band luminosity func¬ 
tion para meters and the direct ly measured parameters 
in g from iLovedav et all 1 2012 ) is good at z < 0.5; and 
(c) while no A GN with z > 0.5 are considered in this 
work, the lFaber et al.l ( 2007 1 study provides a more ro¬ 
bust and convenient way to parametrize the luminosity 
function to redshi fts bey ond z = 0.5 than an extrap¬ 
olation of the ILovedav et al.l (l2012f) measurements, en¬ 
abling us to estimate the numbers of background galax¬ 
ies reliably. For the i b and, we adopt a modi fication of 
the evolution derived bv ILovedav et all (12012H . The M* 
follows the parametric solution found by these authors, 
but we fix its value at 2 > 0.6 to the 2 = 0.6 value. 
For (j >*, we find good agreem ent with our galaxy num¬ 
ber counts by using the ILovedav et alJ ( 2012 ) value at 
2 = 0 and adopting t he redsh ift evolution for the r band, 
noting that the ILovedav et al. ( 2012 ) parametrized solu¬ 
tion found no redshift evolution in this band . The M* 
and <j>* evolution in i found bv ILovedav et alJ (12012) and 
extrapolated to 2 > 0.5 is shown in Fig. [9l For a, the 
values from these various studies are (-1.0, -1.3, -1.05, 
-1.12) for (u, B, r, i). We extrapolate the highest redshift 
values from each study to all higher redshifts for our cal¬ 
culations, although at the depths of our survey, we are 
generally not sensitive to galaxies with 2 > 1.5. 

4. CORRELATION FUNCTIONS 


Filter mii m 

Mi lm - M* 1 

Sample 1 

u 

22.82 

-1.1 

9 

21.95 

-0.8 

r 

22.14 

0.3 

i 

22.06 

1.6 

Sample 2 

u 

23.31 

-0.6 

9 

22.92 

0.2 

r 

23.01 

1.1 

i 

22.81 

2.3 


1 For z = 0.5 


culated at scales that span these regions, even with a 
large simulated set of random points. To verify this, we 
performed a series of simulations with 5000 random data 
points created with regions identical to those in the real 
data. With the random data, we expect zero correlation 
but find that there is significant structure in the corre¬ 
lation function when the masks are applied, and so we 
calculate the correlation functions in the real data on 
scales that minimize the effect of the masks introducing 
structure into the correlation functions that is unrelated 
to the galaxy distributions. 

To perform the correlation function calculations to a 
uniform depth, we confine galaxy magnitudes to the limit 
defined by the shallowest field in the sample. These 
depths are listed in Table [6] for the two cases we con¬ 
sider: Sample 1 is the sample of all quasar fields with the 
HS0624+6901 g band excluded due to its shallow depth 
{g = 21.79) and Sample 2 is a sample of deeper fields, 
constructed by leaving out other shallow fields in each 
filter to achieve the depth listed in Table [6] The Sample 

1 and 2 limiting magnitudes are also shown as vertical 
dotted and dashed lines in Figured In our calculations, 
we will consider other AGN in the primary sample fields 
with redshifts up to 2 = 0.5. Thus, in Tabled we also 
list the absolute magnitudes relative to M* reached for 

2 = 0.5 galaxies in Samples 1 and 2 in order to demon¬ 
strate the coverage of the galaxy luminosity function at 
the redshifts of the AGN in the overall sample. 


4.1. Galaxy-Galaxy Clustering 

We calculate the angular correlation function of the 
galaxy-galaxy sample, 

u(d) tine = A sg 9 ^ (2) 

using the estimator of lLandv &; Szalavl (I1993H : 


^(^)obs — 




2jVgr (0) 

iv rr (0) 


N ir (9) 


(3) 


To examine the galaxy autocorrelation function and 
the galaxy-quasar cross-correlation signal in the 90Prime 
fields, we use a chip-by-chip approach to the calculations. 
Due to the large interchip gaps in the 90Prime mosaic, 
the masked region in chip 3, and the often masked am¬ 
plifier 4 on chip 2 (see Fig. [2J, the correlation functions 
described below contain significant structure when cal¬ 


where N gg and N TT are the pairwise distances between 
all galaxies in the sample and all random positions in 
the survey area at a given angular separation, 9 , normal¬ 
ized by the total numbers of galaxy-galaxy and random- 
random pairwise distances. The N gI term is the normal¬ 
ized number of pairwise distances between all galaxies 
and random points in each 9 bin. 
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The random points are generated by Monte Carlo sim¬ 
ulations which place galaxies at random positions within 
the areas covered by the 90Prime images. These catalogs 
of random positions are masked in the same way as the 
galaxy catalogs. We generated 10 realizations of each 
field in each filter, with ~40,000 points in each realiza¬ 
tion. The Poisson error in each bin is given by 


A w(0) 


1 + 

a/W 


( 4 ) 


This estimator of the angular correlation function must 
be corrected for the integral constraint, which requires 
the function to reach zero at the survey edges, so that 


w(0)true — w(d)obs + C (5) 


where 



uj(0)d£ldCl 


EA^r(fl)4^ (1 - 7) 

E^rr (0) 


( 6 ) 


(|Groth fe Peebleslll977l : IlnfantelH99411 . 

In each filter, we solve for A gg and 7 independently 
using a least squares technique and derive the uncertain¬ 
ties in the parameters by finding the values of each that 
yield xLn + 1 - 


4.2. Galaxy-Quasar Clustering 

For galaxy-quasar clustering, the autocorrelation func¬ 
tion does not app l y, so w e follow ICroom fc Shanks! (|1999d 
and lBrown et al.l (1200111 

w(0)true = w(0) o bs + C = ^ | — 1 + C (7) 

where N gq and N qi are the numbers of galaxy-quasar and 
quasar-random pairs at each value of 0 , normalized by 
the numbers of galaxies and random points respectively 
and C is defined as in Equ. [G] 

Here we solve for A gq for a fixed value of 7 and again 
derive the uncertainties in the parameters by finding the 
values of each that yield Xmin + 1 . 


4.2.1. Correlation Length 

The galaxy-galaxy and the galaxy-quasar angular cor¬ 
relation functions can be used to derive the correlation 
length, r 0 , of the spatial correlation function 

(i+*r (3+e) (8) 


for galaxy-galaxy and ga laxy-q uasar clustering using 
Limber’s equation (|Limbeilll953H . which relates rg to the 
amplitude of the angular correlation function given in 
Equation [2] 


- _ , -yf<r F ( Z ) D 0 1 1 \z)N x {z)N g {z)g{z)dz 
f 0 °° N x (z)N g (z)dz 


where 


k = \/{Tr) 



( 9 ) 

( 10 ) 



Fig. 10.— Number of galaxies in u , g, r, i frames in all 90Prime 
fields from integration of a Schechter luminosity function with pa¬ 
rameters as shown in Figure [9] integrated to the limiting magni¬ 
tudes listed in column 3 of Tabic [5] Solid lines show Sample 1 and 
dashed lines represent Sample 2. 


fixed in physical or comoving coordinates, or e = 7 — 1 
for clustering according to linear theory. Finally, 

g(z) = —[(1 + z) 3 (l + Q m z + Ha ((1 + z) 2 —l) 2 ] ( 11 ) 
c 

and N x corresponds to either N g or N q depending on 
whether we are considering galaxy-galaxy or galaxy- 
quasar clustering. 

The redshift distribution of the galaxies in our 90Prime 
Samples 1 and 2, calculated from the Schechter luminos¬ 
ity function parameters shown in Figure [9] and the mag¬ 
nitude limits determined by the values listed in Table [G] 
in each of the four bands is shown in Figure [TUI 


4.2.2. Correlation Amplitude 

We also calculate the amplitude of 
quasar correlation fu nction , B ?e usin; 
nique d escribed bv lEllingson et all (Il99ll 
(1200 ID, IMcLure fc Dunlon? (|2001bl i. and 
Ra mos Almeida et al.l ( 20131 . 


the galaxy- 
; t he tech- 
, lFinn et all 


recently by 
We first begin with the 
angular correlation function, which describes the clus¬ 
tering of galaxies about some particular object using the 
following relation: 


n(6)dQ = N g [l + w(0)]d£l. (12) 

Here, n(9)dfl represents the number of galaxies within a 
solid angle dH at an angular distance 6 from the object. 
The term N g represents the average galaxy surface den¬ 
sity, and w(6) = A gq 6 1-7 represents the angular cross¬ 
correlation function discussed above. To calculate B g9 , 
we first determine the amplitude, A gq , from the integral 
of the angular correlation function within a fixed angular 
radius, 9, 


A 


gq ~ 


— - 


(13) 


and F(z) describes the redshift evolution of the correla- where N t represents the quasar field galaxy counts, and 

tion function in Equ. [ 8 ] with e = 0 or 7 — 3 for clustering Nt represents the background galaxy counts both within 
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angular distance 6 of the quasar. For each quasar field, 
we count the galaxies within a specified radius around the 
quasar, down to the magnitude limit specified by the am¬ 
plifiers on the chip, mii m . To construct the control fields, 
we use all other regions of the mosaic of the reference 
quasar as well as any non-quasar regions of the images 
of our other quasar fields. We place grids of overlap¬ 
ping control fields with the same radius as our reference 
quasar in each of these regions and count the galaxies 
brighter than mii m . The galaxy counts in these regions 
may be affected by other galaxy clusters and large scale 
structures in the quasar fields. We discuss these in de¬ 
tail for each quasar field in Appendix [X] To mitigate this 
problem we take a median of all the counts from these 
regions and use this as our value for N b . 

The angular covariance amplitude, A gq , quantifies any 
excess in galaxy counts within the quasar field as com¬ 
pared with the background counts of the control field. 
To de-project the angular cross-correlation function into 
its spatial equivalent, we use: 

n(r)dV = p g [l + t,(r)\dV, (14) 

where n(r) and p g represent the total and expected av¬ 
erage galaxy counts within volume dV, respectively, and 
the spatial cross-correlation function, £(r), is defined as 
in Equ. [H with amplitude B gq \ 


&r) = B gq r~ J , (15) 

where 7 = 1.77 ( Groth fc Peebles 1977 ). B gq can be 
obtained using the following equation relating the spatial 
clustering amplitude to th e angular clustering amplitude 
(iLongair fc Seldneilll979l f : 


B 


gq — 


NgA gq 

4*(miim? z'jl'y 


D 


7-3 


1 + 2 


(16) 


Here, D is the angular diameter distance to the quasar, 
and 4>(z) is the galaxy luminosity function integrated to 
the limiting absolute magnitude of the data at the quasar 
redsliift. The quantity / T is a constant of integration 

equivalent to r( || 7 1 ) /2) • 

The errors in A gq and B gq can be calculated using the 
following equations: 


^■A gq _ A B gq 

A Ft 

^gq D gq 


1 

N t -N b 


[(N t -N b ) + 1.3 2 N b )K 


(lYee fc Lopez-Cruzl I1999H . We normalize the galaxy 
counts using the luminosity function calculated with the 
redshift-dependent parameters shown in Fig. [9] 

To calculate the predicted average galaxy background 
counts, N g , we integrate the luminosity function using 
the same redshift-dependent parameters: 




N g = 


4>(to, z ) dmdz 


' 0 J —00 


(18) 


to the relevant luminosity dictated by the limiting mag¬ 
nitude of the data at each redshift. In calculating each 
value of Mi; m , we apply a median K-correction for galax¬ 
ies in each filter at each redshift calculate d from the 
KCORRECT code dBlanton fc Roweisl 120071) and using 
the galaxy templates of lKinnev et al.1 (1199611 . 


In our data sample, we adjust the galaxy magnitudes 
for extinction using galaxy extinction E(B-V) values ac¬ 
quired from the NASA/IPAC Infrared Science Archive 
for each qua sar, and assuming Ry = 3.1 and the redden¬ 
ing curve of lO’Donnelll •: 1 99 1 1 . In our calculations, we 
used a radius of 0.5 h~ l Mpc around the sample quasars 
and for the control fields derived from the other 90Prime 
images of comparable depth in the same fil ter, and we fix 
the parameter q to be 1 . 77. We follow iFinn et, all (120011 1 
and IMcLure fc Dunlod (l2001bl l in using only galaxies in 
the range m*-l - m *+2 in order to eliminate many galax¬ 
ies in the foreground or background of the quasar of in¬ 
terest. 


5. RESULTS 


5.1. Comparisons with other work 
5.1.1. Correlation Amplitude 


Results of the calculations of B gq are listed in Table [7] 
and plotted versus quasar redshift in Figure [Til along 
with values for some of our fields present in the liter¬ 
ature. We do not include a value for PG0844+349, as 
its low redshift, 2 = 0.064, results in an angular radius 
corresponding to 0.5 h^ 1 Mpc that is larger than our 
contiguous chip areas. We do include this object in our 
discussion of the angular correlation function below. We 
also omit the i band data of Q2141+175 from both this 
calculation and the correlation function calculations as 
the quasar was placed in amp 4, which was masked due 
to noise. 

The limiting magnitudes of the galaxies included in 
these calculations are those listed in Table [5] with the 
exception of PG0953+414 u and i bands (23.7 and 23.2, 
respectively), PG1404+226 i (22.5), and PG1612+261 
g (23.5), as there were no galaxies found in the control 
fields in these bands down to the quoted magnitude lim¬ 
its. This change in the magnitude limits had no effect on 
N t . the quasar field galaxy counts, but simply increased 
N b , the median number of background galaxies found in 
the control fields to a number greater than zero, result¬ 
ing in a finite value of B gq . No significant differences 
in the B gq values were found when imposing a uniform 
magnitude limit on all the fields, so we opt for using the 
individual values for each field. All B gq values are quoted 
in units of Mpc 7 , where 7 = 1.77 unless otherwise noted. 
The mean values of B gq we find in ( u , g , r, i) are (21+20, 
161+96, 58+100, and 445+271). 

Gi ven the well-known morphology-density relation 
(jDressleri Il980f ) and the fact that radio loud quasars 
are more commonly found in massive early-type hosts 
(|Hamilton et all [2002t iBest et all 120051) . it may be ex¬ 
pected that galaxies cluster more strongly around radio 
loud quasars than around radio quiet objects. Opti¬ 
cal and X-ray selected AGN show no tendency to re¬ 
sidein galBxy_grou 2 _enyironments (Coldw ell fc Lam bas 
1 2 0031 iSilverman et all I2009H but iWold et al.l (1200+ 


find that radio loud quasars do prefer group or 
Abell class 0 environments. Some other studies have 
found this distinction between the environments of 
radio loud and radio quiet quasars (lEllingson et al. 


1 991 : IBest et al.l 120051 iShen et al.1 I2009T Hickox et al 
2009 : iRamos Almeida et ah 2013T) while others have not 


dFisher et al.1 1996^McLure fc Dunlor 1200 Hi [Finn et al.1 


120011 : iWold et al.ll2001l : lColdwell et al. 200211 . As we have 
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only 2 radio loud quasars in our sample, we have little 
leverage to address this question. However, excluding 
the two radio loud objects in our sample, PG1545+210 
and Q2141+175, we find a systematic decrease in all four 
filters: < B g9 >= (9±18, 144±114, -39±56, 295±260). 
Removing the two additional u fields (PG1444+407 and 
HS0624+690) and one g field (HS0624+690) that do not 
have limiting magnitudes that reach M*+2 at the red- 
shifts of the AGN does not affect these means signif¬ 
icantly, bringing them to 32±14 and 202±113 respec¬ 
tively. The error bars reflect the standard deviation in 
the measurements among all fields in each filter. The me¬ 
dian redshift of our sample AGN is 0.18 and the median 
M r is -23.7. These modified B gg averages are all con¬ 
sistent with th e galaxy-galaxy correla tion amplitude, 20 
(ft -1 Mpc) 1 ' 77 (iDavis fc Peeblesl 1 1 9831 ) . although within 
the uncertainties, the average values in g and i are also 
consistent with richer environments, as discussed in more 
d etai l for individual fields below. 

iFinn et ahl (120011) used images of quasar fields taken 
with the F606W, F675W, and F702W filters and the 2.5' 
field of view of the Wide Field Planetary Camera 2 on 
the Hubble Space Telescope and, scaling to our assumed 
value of H 0 , find <B gg >= 104±22, -33±27, and 133±44 
in those respective filters for radio qui et quasars at me- 
dian redshifts of 0.17, 0.42 and 0.4. IMcLure &; Dunlopl 
(|2001bf) find <B gg >= 326 ± 94 for their sample of radio 
quiet objects witfi median redshift of 0.16, also derived 
from F606W and F675W /LS T/WFPC2 im a ges. Their 
dataset is a superset of that of [Fisher et al.l (119961 ) who 
found <B gg >= 132±37. These HST studies were limited 
to radii of ~200 kpc, while ground-based studies typi¬ 
cally extend to the same 0.5 ft -1 Mpc scale as ours does. 
The g round-based V, R, and I photometry of lWold et all 
( 2001 1 targets AGN at higher redshift than our AGN 
sample, z me d = 0.714, and they find <B gg >= 210 ± 82 
for radio quiet quasars. All of these studies show av¬ 
erage B gg values that are larger than B gg at the 2 — 3a 
level. Ou r overall resu l ts are mor e in line with the earlie r 
studies of lSmith et al.1 (11995| ) and lEllinerson et ahl (|1991h . 
both with z me d ~ 0.4 who found average B gg values con¬ 
sistent with B gg for radio quiet objects. 

Considering the individual values, we find that there 
can be a large range of B gg values obtained for the 
same field depending on the filter used. Typically, the 
values of B gg in the i band tend to be larger than 
those of the other bands and those calculated in the 
u band are nearly all consistent with B gg . As above, 
for the r band, we compare our indi vidual resu lts with 
thos e found in the literature, from IFinn et ahl (|2001f ) 
and IMcLure fe Dunlopl (1200 lH) . and for PG1545+210, 
we compare with lYee fc Greenl (|1984T l. who calculated 
B gg from r band images taken with the SIT-vidicon cam¬ 
era and the 1.52 m Palomar telescope. We scaled these 
literature B gg values for our value of Ho but otherwise 
do not corre ct for different cosmologie s used in these 
other works (fMcDure fe Dunlopl (l2001bh : Hq = 50 km 
s -1 Mpc -1 : [Finn et ahl (120011) Hq = 100 km s -1 Mpc -1 , 
both use = 1, Qa = 0). In most cases our results are 
consistent with these previous studies. We give a more 
detailed discussion of these comparisons for each field in 
A ppendix lAl 

lYee & Ellingsonl (119931 ) considered B gg > 500 to 


TABLE 7 

B gg Measurements 


QSO 

Bgg (Mpc 7 ) 1 

Filter 

Ref. 2 

MRK586 

-3±6 

u 

1 


-36F137 

g 

1 


0 3 

r 

1 


93F161 

F606W 

2 


13±82 

F606W 

3 

HS0624+6907 

-66±52 

u 

1 


-319F279 

g 

1 


-144T127 

r 

1 


-126F151 

i 

1 

PG0923+201 

-76±25 

u 

1 


-114±133 

g 

1 


-377F130 

r 

1 


-529±291 

i 

1 


500F269 

F675W 

2 


119±81 

F606W 

3 

PG0953+414 

30±28 

u 

1 


207F101 

g 

1 


-36F191 

r 

1 


629F519 

i 

1 


730F297 

F675W 

2 


159±81 

F606W 

3 

PG1116+215 

68±41 

u 

1 


-41±57 

g 

1 


186F134 

r 

1 


863F461 

i 

1 


321±235 

F606W 

2 


185±81 

F606W 

3 

PG1307+085 

10±27 

u 

1 


259F120 

g 

1 


125F306 

r 

1 


1461F778 

i 

1 


112±176 

F606W 

2 


170±82 

F606W 

3 

PG1404+226 

85±85 

u 

1 


942±522 

g 

1 


143F759 

r 

1 


707F854 

i 

1 

PG1444+407 

-27±19 

u 

1 


46F119 

g 

1 


-135F161 

r 

1 


292±334 

i 

1 


57F130 

F606W 

2 


49±79 

F606W 

3 

PG1545+210 

133±38 

u 

1 


104±120 

g 

1 


0 3 

r 

1 


1643F404 

i 

1 


206F199 

F606W 

2 


242±79 

F606W 

3 


236F120 

r 

4 

PG1612+261 

56±27 

u 

1 


354±237 

g 

1 


-116F397 

r 

1 


-935F727 

i 

1 

Q2141+175 

370F221 

g 

1 


991F355 

r 

1 


112±169 

F675W 

2 


3 Here we quote results for 7 — 1.77 and a radius of 
0.5 /i - 1 Mpc; Literature values have been scaled by 
(h ~ 1 ) 1 ' 77 , but are otherwise uncorrected for different 


cosmologies. 

2 References: (l)this work; (2 lMcLure fc Dunlopl 

Il2001bll : O lFinn et al.l ll2Mill : )4 1Yee fc Greenl 

IT984D 

3 N t = N b 

be a rich cluster environment, and clusters with 
Abell classes 0, 1 and 2 show values of ~35 0, 

qo 650, an d ~950 respectively , accord ing to lYates et~ahl 
(|1989f) . IMcLure fe Dunlopf (l2001bl) suggest a classi¬ 
fication based on that of lYee fc Lopez-Cruzl dl999l) . 
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Fig. 11.— B gq versus quasar redshift for all filters. Points have 
been offset slightly in z for clarity. Black points show values from 
the literature listed in Table (71 Dott ed lines, bottom to top, m ark 
the B gq values for Abell classes 0-4 dYee Sz L6pez-CruzH1999T) and 
the dashed line marks the galaxy-galaxy correlation value from 
I Davis &; Peebles! (119830 . 

where B 59 =(300,700,1100,1500,1900,2300) corresponds 
to Abell class (0,1,2,3,4,5) respectively. Given the large 
statistical uncertainties in the B gq values, we find re¬ 
sults marginally greater than 500 in only a few cases, for 
PG1307+085 and PG1545+210 in the i band, 1461±778 
and 1643±404, respectively, for PG1404+226 in the g 
band, 942±522, and for Q2141+175 in the r band, 
991±355. The radio loud AGN PG1545+210 is known 
to re side in an Abell class 1 environment (lOemler et al.1 
Il972h . though our i band B gq value is consistent with 
Abell class 3. The B gq values for this field in g and r 
are consistent with B gg , and in u we find only 133±38, 
which is significantly larger than B sg , but far smaller 
than the i band value. PG1307+085 has a foreground 
galaxy cluster in its field at a close projected separation 
from the quasar and with a similar redshift, z c i us = 0.14 
versus z q = 0.155; and PG1404+226 has a background 
cluster within the 13. 1' radius corresponding to 0.5 /i -1 
Mpc. Q2141+175 is the only other radio loud AGN in 
our sample and while it has no known cluster at its red¬ 
shift or in projection, it is a less well-studied field than 
many of the others. Thus the B gq for these other fields 
may be influenced by foreground, background, or even 
proximate galaxy clusters, but only in select bands. 

There are a handful of other fields that show a sig¬ 
nificant, positive signal in B gq in one or more bands, 
and these are discussed in Appendix [A] One field, 
PG0923+201, shows negative values of B gq in all bands, 
with two (u,r) at > 2cr. As we discuss in Appendix lAl 
this well-studied field has a compact group consistent 
with the redshift of the AGN and with a center only 23 
h~ x kpc from the position of the AGN, making this re¬ 
sult an intriguing one. Three other fields, HS0624+690, 
MRK586, and PG1444+407 show B gq values in all filters 
that are either consistent with or significantly less than 
B gg . HS0624+690 is the only quasar in our sample, and 
with z = 0.370, our limiting magnitudes in u and g allow 
us to reach only ~ M*+0.7 and ~ M* 3-0.1, respectively. 


TABLE 8 

Galaxy-Galaxy Clustering 
Clustering in 90Prime QSO fields 1 


Filter 

7 

^gg 

Xu 

u 2 

1.513=0.02 

0.133=0.05 

2.9 

9 2 

1.50±0.03 

0.263=0.08 

10.9 

r 

1.59±0.03 

0.183=0.04 

12.5 

i 

1.533=0.02 

0.183=0.03 

27.1 


1 0.1667'- 10.1667' 

2 Excludes chips with one or more ampli- 
fiers masked 


Two other AGN in our sample, PG1116+215 and 
PG1612+261 lie near galaxy clusters, though their red- 
shifts and position on the sky place them near the pe¬ 
ripheries of those clusters. The B gq values for these 
fields give no indication of a cluster association: for 
PG1116+215, all the B gq values are consistent with B gg , 
except the anomalously high value we find for the i band, 
8633=461. For PG1612+261, the only filter showing a B gq 
marginally larger than B gg is the g band, for which we 
find 354T237, Abell class 0. 

5.1.2. Correlation Functions 

We show the results for the normalization, A gg , and 
slope, 7 of the galaxy autocorrelation function in Fig¬ 
ure [12] and list the results in Table [ 8 ] These calculations 
were performed for Sample 2. For the u and g filters, 
we use only chips 1 and 4 as chips 2 and 3 added extra 
structure due to masking. The r and i filters typically 
are less affected by this. The use of all chips for u and 
g gives similar results, but a poorer fit to the form of 
the autocorrelation function. We avoid separations less 
than 10 " to mitigate against blending, and because we 
are restricted to single chips, we calculate the correlation 
function to a maximum separation of 10'. The errors on 
the parameters are derived from the variance in that pa¬ 
rameter from 100 jackknife samples of the dataset, where 
the sampling with replacement was done one field/filter 
at a time. 

In our calculations of the galaxy-quasar cross¬ 
correlations, we consider quasars in 90Prime fields other 
than the target quasars themselves, as any clustering 
may also be present for these objects as well. As dis¬ 
cussed above, the galaxy counts in our images are very 
sparse at z > 1 and so galaxies associated with clus¬ 
ters, groups, or quasars at these redshifts are generally 
not of concern to us here. AGN with z < 0.5 found in 
the 90Prime frames are listed in Table |9] and the result¬ 
ing quasar redshift distributions in each filter are shown 
in Figure [13] The absolute magnitudes relative to M* 
reached for z = 0.5 galaxies in Samples 1 and 2 are listed 
in Table 0 

We list the results for the normalization A gq of the 
galaxy-quasar cross-correlation function in Table [TO] for 
both Samples 1 and 2, where the latter correspond to 
deeper magnitude limits and fainter galaxy samples. For 
these solutions, we hold the slope 7 to be constant at 1 . 8 , 
the commonly accepted value and 1.5, corresponding to 
the slope we found from the 90Prime data for the galaxy- 
galaxy clustering on 10' scales. This difference in slope 
changes the normalization at only the ~ 1 <t level. The 
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Fig. 12.— Galaxy-galaxy angular correlation function in u, g, r, i 
filters for 90Prime fields. 



Redshift 

Fig. 13.— Redshift distributions of all quasars in the 90Prime 
fields 


TABLE 9 

Other z < 0.5 Quasars in 90Prime Fields 


QSO 

z 

Filters 

PG0844+349 

SDSS J084731.77+351416.4 

0.237 

u , r 

2MASX J08503620+3455231 

0.144 

u , r 

PG0923+201 

SUSS J 092507.72+203540.9 

0.472 

9,i,r 

SDSS J092525.16+202139.0 

0.460 

u , g , r 

SDSS J092536.08+201649.5 

0.228 

9 

PG1307+085 

Q87GBBWE91 1307+0843 


r , i 

SDSS J131155.76+085340.9 

0.469 

u, g, r, i 

PG1444+407 

SDSS J144618.08+412003.0 

0.268 

u,g,r,i 

PG1545+210 

SUSS .1154749.70+205056.5 

0.265 

u, g , r 

SDSS ,1154750.71+210351.1 

0.296 

u , g , r 

SDSS J155014.81+212431.5 

0.479 

r 

SDSS J155046.30+205803.1 

0.401 

u , Q , r\ i 

2MASS J15505930+2128088 

0.372 

u, r 

PG1612+261 

HB89 1612+266 

0.395 

u, g, r, i 


r esults for both slopes are plotted in Figure fl4l 

iCroorn fc Shanksl (|1999T) found no galaxy-quasar cor¬ 
relation in their sample of ~150 optically and X-ray se¬ 
lected quasars with z = 1 — 1.5 and bj < 23 galaxies in 
five deep A AT plate fields. In fact, their data showed a 
weak anti-correlation between these two populations. At 
the the lower redshifts of our sample, we see significant 
correlations only in the u band. For Sample 1, there is 
also a significant positive correlation on the r band, but 
this does not hold with the inclusion of fainter galax¬ 
ies in Sample 2. In fact in all bands except the i band, 
the correlation is significantly smaller for Sample 2, and 
we also see weak anti-correlations in several cases. In i, 
the normalizations are all consistent with zero for both 
Samples 1 and 2. 

In Table ITOl we also tabulate the correlation length for 
galaxy-quasar clustering, ro, found from Limber’s equa- 


TABLE 10 

Galaxy-Quasar Clustering in 90Prime 
QSO fields 1 


Filter 

^gq 


Xv 

r 0 h” 1 Mpc 2 


Sample 1, 7 = 

1.8 

u 

0.49Lri4 


0.39 

8.07 

9 

0-23±g;it 


0.69 

5.36 

r 

0 21 + 0 - 08 
U.Z1_Q 07 


0.70 

6.11 

i 

-0.03+0.07 


1.66 



Sample 

1 , 7 = 

1.5 

u 

u - Oi -0.19 


0.38 

16.9 

9 

0 34 + 0 - 22 
U-c^_o.i9 


0.65 

11.64 

r 

0.28+0.10 


0.65 

12.3 

i 

-0.0027+i° 


1.68 



Sample 

2 , 

7 = 

1.8 

u 

0-40t£u 


0.61 

7.76 

9 

-0.04+0.09 


1.03 


r 

-0.05+0.05 


1.57 


i 

0.004+0.05 


1.44 

0.90 


Sample 

2 , 

7 = 

1.5 

u 

o.sitrus 


0.57 

16.2 

9 

-0-03+1? 


1.04 


r 

-fi Ofi ” 1 ” 0,07 
-U.UO_o 06 


1.58 


i 

0.02+0.07 


1.44 

2.95 


1 0.1667'- 10.1667' 

2 e = 0 




fl(arcmin) fl(arcmin) 

Fig. 14. — Quasar-galaxy angular correlation functions in 
u,g,r,i filters for 90Prime fields for 7 = 1.8 and 7 = 1.5. Points 
are offset in x and g fit curve is offset in y by + 0.02 for clarity. 
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Fig. 15.— Comparison of the two methods for calculating the an- 
gular covariance amplitude, A gq from Equation [7] and Equation ll31 
both with 7 = 1.8. 

tion, Equation [51 and the galaxy and quasar redshift dis¬ 
tributions plotted in Figures [lU] and [13] This is undefined 
for negative values of A qq , so no correlation length values 
are reported in Table [10] in these cases. The values we 
find range from 0.9 ft ” 1 Mpc for Sample 2 in the i band 
with 7 = 1.8 to ~16 ft ” 1 Mpc for both Samples 1 and 
2 in the u band with 7 = 1.5. These results are gener¬ 
ally within the range found for 0 . 7 < z < 1.4 quasars in 
the DEEP2 survey bv lCoil et al.1 (1200711 . 0.1 ft ” 1 < r 0 < 
10 ft ” 1 Mpc. iBrown et al.1 (1200111 studied a population 
of galaxy fields around a sample of moderate redshift 
quasars with z = 0.2 — 0.7 and also find a weak corre¬ 
lation, though a stronger one for red galaxies (ro ~ 7 
ft ” 1 Mpc, 7 ~ 1.9) than for blue galaxies (r 0 < 4 ft ” 1 
Mpc). These authors provide a summary of results to 
that time in their Table 1. The value of the correlation 
length ranges from ~5-8 ft ” 1 Mpc for samples of radio 
quiet quasars with similar redshift but shallower depth 
than our current 90Prime study to ~12-17 ft ” 1 Mpc for 
fields around radio loud quasars, comparable to our re¬ 
sult for the u band in the 90Prime sample. Below where 
we discuss trends of galaxy-quasar clustering with quasar 
properties, we look at the question of whether the two ra¬ 
dio loud objects in our sample show significantly greater 
galaxy clustering than the radio quiet objects. 

In our two methods of investigating galaxy-quasar cor¬ 
relations, outlined in Sections 14 21 and [4. 2.21 we calculate 
the angular covariance amplitude, A gq in two different 
ways, one involving a comparison with a set of random 
galaxy positions and one using the expectation from an 
integral of the galaxy luminosity function in the relevant 
band. In Figure IT5l we compare the values obtained from 
these two methods for 7 = 1.8. We note that the uncer¬ 
tainties in the fitted uj(9) for a single field/filter can be 
large, and the offset for any one field/filter can be large, 
ranging from - 0.6 to +1.8, but the mean offset in these 
values is 0.22 with a standard deviation of 0.55. As has 
been noted in previous work, clustering statistics like B gg 
can be sensitive to the methodology used to calculate it, 
the choice of luminosity function and control fields. We 


TABLE 11 

Sample Quasar Properties 


QSO 

M r 

log(M B H) 

Refn. 1 

[O III] Luminosity 
erg/s xlO 42 

Refn. 2 

MRK586 

-23.4 

8.34 

1 

3.50 

1 

HS0624+690 

-26.0 

9.69 

2 



PG0844+349 

-22.4 

7.96 

3 



PG0923+201 

-23.7 

8.00 

4 

0.99 

2 

PG0953+414 

-24.7 

8.44 

3 

1.90 

3 

PG1116+215 

-24.6 

8.52 

4 

2.60 

3 

PG1307+085 

-23.2 

8.64 

3 

3.17 

2 

PG1404+226 

-21.8 

6.88 

4 

0.38 

2 

PG1444+407 

-24.1 

8.28 

4 

2.22 3 

2 

PG1545+210 

-24.4 

9.31 

4 

4.93 

2 

PG1612+261 

-22.6 

8.05 

4 

6.68 4 

2 

Q2141+175 

-23.6 

8.74 

1 

2.80 

1 


4 References: ( l lMr.Lure_fe^Dunlorj -d2nniafI: f2 1Labita et all (1200011 : 

M l Peterson et all 110011) : ( 4 1 Vestergaard fc Petersonl (1200611 
2 References: (T lHo fc Kiml (1200011 : (2) This study; 13 1Shang et all (1200711 
^ IShang et al.1 (1200711 find 0.22 

4 Sum of two components FWHM=420 km s —x and 1000 km s —1 with fluxes 4.75 
and 1.93 respectively 

therefore concentrate less on the absolute values of this 
parameter than on investigating any trends in the data, 
discussed in the following section. 

5.2. Trends in Clustering Parameters 

5.2.1. Correlation Amplitude 

We investigate trends of galaxy-quasar clustering 
with quasar absolute r magnitude, black hole mass, 
and [O I llj line luminosity, as an indicator of AGN 
activity dKauffmann et al.l 120041) . These properties 
are listed in Table .IT The values of M r are de¬ 
rived from the r listed in the SDSS DR9 catalog, 
with Galactic ext incti on calcul ated us ing the redden¬ 
ing values from iSchlaflv fc Finkbeinerl ( 2011 1 and the 
lO PonneTil (I1994T) reddening law. We also applied a 
K-corre ction, calculated from the QSO composite spec¬ 
trum of IVanden Berk et al.l (1200111 and the SDSS filter 
response curves. We collected black hole mass and [O III] 

A5007 line luminosity measurements from the literature 
and list them in Table [TTl Where the [O III] line lu¬ 
minosities do not exist in the literature, we downloaded 
the SDSS DR10 spectra for these objects and fit the H (3 
[O III] A4959,5007 complex using the IRAF task specf it 
and list the resulting value in Tabic fTTl 

Plots of these quasar properties with P> gq parameter 
are shown in Figures I16H181 but in no case do we find a 
significant trend of this clustering measure with any of 
the quasar properties. 

5.2.2. Correlation Functions 

We also investigate trends of the galaxy-quasar angular 
correlation function with quasar luminosity, black hole 
mass, and [O III] line luminosity by using the relative bias 
with respect to galaxies, b gq = (uigq/wgg) 1 / 2 . Overall, 
we find b gq = (3.08 ± 0.51,1.49 ± 0.53,1.48 ± 0.29, 0.95 ± 

1.49) for (u, g, r, i) for cross- and autocorrelation function 
solutions with fixed 7 = 1.5. The bias is only significantly 
greater than unity for the u band and decreases with 
increasing wavelength, perhaps reflecting a trend seen 
for SDSS Sy 2 galaxies to lie in bluer environments than 

















































































18 



-26 -25 -24 -23 -22 

M r 


Fig. 16. — B gq in each filter versus quasar absolute r magni¬ 
tude. Solid lines and shaded regions mark the mean and dispersion 
among all the values in a given filter, dotted line s show t he expected 
corre lation amplitude for Abell classes 0-4 dYec & Lope z-Cruzi 
1999), an d the dashed line mark s the galaxy-galaxy correlation 
value from I Davis &; Peebles! (j 19831) . Black points show values from 
the literature listed in Table [7] 



Fig. 17.— As in Figure fT6l but for B gq in each filter versus quasar 
black hole mass 

normal galaxies (|Coldwell et al.ll2014f ) . This runs counter 
to a trend in the mean B sg values quoted above, for which 
the u band values are lowest and those in i are the largest. 
This is driven by the large i band B gq values found for 
PG1545+201 and PG1307+085, the former being a radio 
loud quasar. That trend is less apparent for the radio 
quiet objects alone, and without PG1307+085 as well, 
the average i band B gq value drops to 129±238. 

The values of ui for each field and filter individually 
are listed in Table m For each field/filter where both 
the cross- and autocorrelation values are positive in the 
innermost bin, to < 1.167', we find no correlations in 
bgq with any of the quasar properties discussed above. 
Most fields for which both m gq and uj gg within 1.167 / 



bom] 5007 10 « ergs/s 

Fig. 18.— As in Figure fl6l but for 3 gq in each filter versus quasar 
[O III] luminosity 

are positive show b gq > 1, but only a handful at > 2a: 
PG0844+349 and PG1612+261 in u, HS0624+690 in g, 
PG0923+201 and PG1116+215 in r, and PG1404+226 
in i. For PG0844+349 and PG1612+261 the filters with 
b gq > 1 are the only filters with u> gq > 0, for which b gq 
is defined. That the u filters would show a large posi¬ 
tive relative bias while all others show a quasar-galaxy 
anticorrelation may be indicative of bluer environments, 
although this does not hold for the HS0624+690 field, for 
which the b gq in g is 1.31 ± 0.094 but the values of ui gq 
are negative in all other filters, including u. A few fields 
show b gq > 1 consistently in all filters (PG0923+201, 
PG1116+215) or in all filters but one: PG1307+085 (i), 
PG1404+226 (r). The bias seen in the u band for the 
radio loud quasar PG1545+210 is negative, and not sig¬ 
nificantly greater than unity in any filter. Similarly, for 
the other radio loud AGN in our sample, Q2141+175, 
only g and r band data are available and neither shows 
a large bias. In g , ui gq < 0, and in r b gq = 0.66 ± 0.75. 

6 . DISCUSSION 

A somewhat contradictory picture emerges when we 
compare the results of the correlation function analysis 
to the calculated B gq statistics. The fields that show 
the largest values of B gq , PG1307+085, PG1545+210, 
PG0953+414 in the i band, PG1404+226 in g , and 
Q2141+175 in r, all show relative biases that are con¬ 
sistent with unity. The PG0923+201 field showed an 
anticlustering B gq signal in u and r, but both filters give 
b g q > 1, with the r band significant at ~ 3cr. We note 
that, in addition to a luminosity function dependence, 
the B g q analysis probes a larger spatial scale than the 
relative bias calculation within 1.167', corresponding to 
~10-40h~ l kpc over the redshift range of the AGN in our 
sample. 

Overall, as noted above in Section 15. 1. 11 the mean val¬ 
ues of the correlation amplitude for 0.5 h~ x Mpc scale 
clustering, while all consistent within the uncertainties 
with galaxy-galaxy clustering, show the largest signal in 
the i band, B 9g =295±260. The correlation function anal¬ 
ysis, on the other hand, probing ~10-40h -1 kpc scales 
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TABLE 12 

Galaxy-Galaxy and Galaxy-Quasar 
Correlations for each 90Prime Field 1 


QSO 

Filter 

u ss 

OJgq 

MRK586 

u 

0.151i0.036 

0.08il.08 


g 

0.015i0.008 

0.33i0.54 


r 

0.072i0.006 

0.07F0.47 

HS0624+6907 

u 

0.071±0.024 

-0.67±0.32 


g 

0.387i0.055 

0.67 2 


r 

0.041i0.008 

-0.50i0.24 


i 

0.082±0.008 

-0.98i0.01 

PG0844+349 

u 

0.215±0.022 

0.95i0.50 


r 

0.006±0.003 

-0.08i0.12 

PG0923+201 

u 

0.150i0.019 

0.19i0.53 


g 

0.194±0.010 

0.57i0.36 


r 

0.223±0.011 

0.95i0.32 


i 

0.168±0.005 

0.54i0.37 

PG0953+414 

u 

0.206±0.012 

0.32i0.47 


g 

0.344±0.020 

1.65i0.93 


r 

0.246±0.008 

-0.36i0.31 


i 

0.270i0.004 

0.19i0.30 

PG1116+215 

u 

0.433±0.028 

1.13i0.86 


g 

0.630i0.039 

1.80il.40 


r 

0.252i0.017 

2.47il.31 


i 

0.127i0.006 

0.83i0.61 

PG1307+085 

u 

0.094i0.023 

0.39i0.62 


g 

0.158±0.015 

0.42i0.53 


r 

0.088±0.007 

0.35i0.27 


i 

0.104i0.005 

0.OliO.19 

PG1404+226 

u 

0.126i0.013 

0.54i0.77 


g 

0.123i0.007 

0.24i0.87 


r 

0.171i0.004 

-0.15i0.28 


i 

0.170i0.005 

1.09i0.63 

PG1444+407 

u 

0.257±0.022 

0.36i0.96 


g 

0.029i0.008 

-0.03i0.30 


r 

0.159i0.007 

-0.35i0.20 


i 

0.115i0.005 

0.10i0.24 

PG1545+210 

u 

-0.029±0.014 

0.33i0.25 


g 

0.257i0.018 

-0.46i0.16 


r 

0.063±0.006 

-0.17i0.14 


i 

0.093i0.005 

0.09i0.21 

PG1612+261 

u 

0.178±0.023 

1.84i0.79 


g 

0.432i0.009 

-0.47i0.18 


r 

0.203i0.006 

-0.47i0.15 


i 

0.199i0.005 

-0.50i0.13 

Q2141+175 

g 

0.172±0.009 

-0.23i0.27 


r 

0.389±0.008 

0.17i0.39 


1 Calcuated within 0.1667 / < 9 < 1.1667* for fixed 7 = 
1.6 

2 No galaxies in this inner bin, so formal error is inde¬ 


terminate. 

shows the largest bias for quasar galaxy clustering in the 
u band, with galaxy clustering around quasars a factor of 
3 larger than the galaxy autocorrelation signal at these 
scales. This suggests that different galaxy types may 
cluster around quasars on different scales. Some ear¬ 
lier studies have seen clustering signals at small scales 
and differences in galaxy populations near quasars. In a 
future paper, we will use the multi-band data to calcu¬ 
late photometric redshifts for the galaxies in our sample, 


and investigate trends in the clustering parameters with 
galaxy types determined with these redshifts. 

The basic trends in our results are consistent with ear¬ 
lier findings that quasars cluster like L* galaxies, and 
indeed like galaxies up to 2 magnitudes fainter than L*. 
We also confirm th at there are no significant trends with 
quasar luminosi ty dFinn et all 120011: iSerber et all 120061 
iCoil et aLll2007ll . though unlike ISerber et al.l (120061) . we 
also find that the lack of luminosity dependence holds 
at small as well as large spatial scales. In our sample of 
broad line AGN, we do not find the trend for AGN with 
lower [O III] luminosities to have la rger c l ustering sign als 
as seen in SDSS narrow line AGN (| Wake et al .11 20 041) . 

Thus, the mechanism driving accretion onto the central 
black holes in these objects is not imprinted on their large 
scale galaxy environments. This may be due to a non¬ 
merger mode in fueling the central black holes or to a 
delay in the onset of nuclear activity from the time of 
the initial galaxy interactions which triggered its fueling. 
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APPENDIX 

NOTES ON INDIVIDUAL OBJECTS 

In this section, we outline notes on each of our 90Prime fields, based on previous studies of the quasar fields and host 
galaxies, searches for nearby objects in the NASA Extragalactic Database (NED), and the results of our calculations 
of the correlation amplitude, B gq . We refer the reader to Table [T] and Figures UTI fT8l 















20 


Mrk 586, z = 0.155 

The field of this radio quiet object has been imaged previou sly in narrow band [O III] by IStockton fe MacKentvl 
( 1987 1 and in the optical (B ,g,i filters) bv lKirhakos et al.l ( 1994 ) using the Palomar 1.5 m (image sizes ~3' x 3') who 
note four galaxies in this field, including a large disk galaxy (NAB 0205+02: 06) at a projected separation of 15.6 /i _1 
kpc. __ 

The quasar host was studied by (|Bahcall et al.lll997h using the HST/ Wide Field/Planetary Camera-2 (WFPC2) and 
the F606W filt er. The host is disk-like and the AGN has a com panion with jet-like structure at a projected distance 
of 23 kpc (|Stockton fc MacKentvl U987t lBahca.ll et al.lll997ll . There is no evidence of tidal features to suggest a 
direct interaction. 

A NED search reveals ~17 galaxy clusters, groups, or candidates, within 50' of the quasar. The cluster with the 
smallest angular separation from the quasar (7.7') is NSCS J020723+024645, which lies at z = 0.37. 

The B gq we calculate from the 90Prime r band image is identically zero because the number of galaxies detected at 
5cr within 0.5 h -1 Mpc of the quasar exactly equ a ls the m edi an number found i n 871 control fields. This is consistent 
with the values found by both IMcLure fc Dunlod (l2001bl) and lFinn et al.1 (120011) from the HST/ WFPC2 images in the 
F606W filter, 93+161 and 13+82, respectively, and with the values calculated from the 90Prime u and g band images, 
-3+6 and -36+137, respectively. We have no 90Prime data of this field in the i band. 

HS0624+ 690, z = 0.370, QSO 

The host of this radio quiet quasar is an elliptical galaxy dFlovd et al.1 120041) . iKirhakos et al.l (Il994f ) identify six 
galaxies in their ground-based optical images within the ~3' x 3' frame. From NED, we find that within 50' of the 
quasar, there is one known galaxy cluster, Abell 0557 (redshift unknown), at 8.03' separation from the quasar. 

All the B gq values we calculate for this field are less than zero, indicating an underdensity of galaxies within 0.5 h~ x 
Mpc of the quasar, however, none of these underdensities is significantly less than B 99 at > 2o\ We conclude that the 
field around this object is consistent with the field. A caveat to this is that HS0624+690 is the the highest redshift 
AGN in our sample, and the limiting magnitudes of the u and g coadded images are not significantly fainter than 
M* at the quasar redshift, M* + 0.7 and M* + 0.1, respectively, so we are not reaching as far down the luminosity 
function for this field as the others in the sample. 


PG0844+349, z = 0.064 

This radio quiet object has a host with bulge+dis k structure, sh owing no bar, but some evidence of morphological 
disturbance from deep H-band imaging (iGuvon et al.1 120061 IVeilleux et al.l I2009D . This object qualifies as having a 
bright interacting compa nion (H band luminosity greater than 1/10 L* and projected distan ce within 24 h~ x kpc 
iMcLeod fe Rieka (1 199411) . though it is not likely currently undergoing a strong interaction since IVeilleux et all (l2009f) 
classif y this objec t as a post-merger object with a single nucleus. Also, from the [O 11 ] measurement of I Wilkes et all 
( 1999 ). iHol ( 2005 ) estimates a star formation rate of the host galaxy of less than 1.6 M Q yr -1 . Within a 50' radius 
around this quasar, NED returns 58 galaxy clusters, groups, or candidates. 

PG0923+201, z = 0.192 

The host of this radio quiet AGN is an elliptical galaxy with no evidence o f a bar or of morpho logical disturbance 
from an interaction with a companion (iDunlop et al.ll 20 (Kl IGuvon et al.! 2006; Veilleux et al.l [2009). Previous authors 
have note the pre sence of two bright galaxies within 16 fe -1 kpc of the quasar and have suggested it is a member 
of a small group (He ckman et al.1 11984: Mc Leod fc R.iekel Il994bt iB ahcall e t al.l 1 19971 ). Indeed, the compact group 
SDSSCGB 30540 (z = 0.19) identified bv iMcConnachie et al.l ( 20091) in the SDSS DR 6 has a redshift consistent with 
that of PG0923+201. While the quasar is not listed as a group member in their catalog, it has a separation of ~23 
(i _1 kpc from the group center. There are ~45 other galaxy clusters or groups within 50' of this quasar. 

Our l+„, values from the r and i band fr ames, -377+130 and -529+291, are inconsistent with those of 
iMcLure fc Dunlonl ( 2001b ) and lFinn et al.1 ( 2001 ) at the 3.5-4<r level. These authors find moderately significant over¬ 
densities in this field while we see underdensities in all filters, of increasing magnitude from u to i. This underdensity 
is surprising if the quasar is in fact associated with SDSSCGB 30540. 


PG0953+414, z = 0.234 

Du nlop .et al] ( 2003;) classify the host galaxy of this radio quiet AGN as an elliptical, while others •: Bahcall et al 


11997 : Guvon et al.ll2006l ) have found it to show disk-like structure but to be t oo fai nt to classify reliably. IKirhakos et al 
( 1994T ) find 22 likely galaxies in their optical images described above in Sec. IA.l| 

There are 54 galaxy groups, clusters or candidates within 50' of this quasar, one of which, GMBCG 
J149.20633+41.27869, has a projected separation of only 1.455', though it lies well behind PG0953+414, at z = 0.362. 
A second cluster, WHL J095714.9+411700, is separated by 4.532' from t he qua s ar an d has a redshift z = 0.446. 

Our i band result for B gq is consistent with the result of IMcLure fc Dunlopl (I2001H ) using HST/ WFPC2 F675W. 

600 — 700, but with large statistical uncertainty, ~ 300 — 400. The g band value for this 


Both values are large, ^ 
field, 207+101, is one of 

90Prime r band image is negative but consistent with B gg , and consistent with the value found by iFinn et ahl (|2001[ ) 
from HST/WFPC2 F606W data within the uncertainties. 


1/4 of the B gq results that are significantly greater than B 9S , while the B ?a _result from the 






















































































21 


PG1116+215, z = 0.177 

PSF probl ems affected the deep H band imaging of the host galaxy of this radio quiet object, with some conflicting 
resul ts. [Guv on et al, (2006! ) find that the host is disk-like and asymmetric with some evidence for a disturbance, while 
iVeilleux et al. ( 20091) find that it shows no spri al arms or bar or m orphological disturbance and classify it as an old 
merger, in agreement with the optical imaging (Bahcall et alJfi~997l ). 

There are ~68 galaxy clusters, groups, or candidates within this field around PG1116+215. Th e two with the smalles t 
separations from the quasar (0.28' and 1.92') are ZwCl 1116.5+2136 and SDSSCGB 16808 (iMcGonnaclrie et al.ll2009 ). 
both with unknown redshifts. However, one cluster found in the SDSS data, MaxBCG J169.85817+21.20845, lies 
at 7.9' separation and has a redshift marginally consistent with PG1116+215, z p h 0 t = 0.170. The relatively large 
projected and redshift separations, ~1 /i” 1 Mpc, and ~1630 km s” 1 , place it toward the periphery of this cluster if it 
is associated with it at all. 

Our results for B gq from the 90Prime r band image, 186+134, is consistent with iMcLure fe Dunlopl (l2001bl ) and 
I Finn et ahl (2001 ) within the large statistical uncertainties. The B gq values found in the u and g filters are also consistent 
with B gg , strengthening the conclusion that PG1116+215 lies in the outer regions or even outside the cluster potential 


of MaxBCG J169.85817+21.20845. 
picture as it is greater than B gg , 


However, the result for the i band, 863+461, 


contradiction to this 


presents a __ 

more consistent with an Abell class 1 environment (IMcLure & DunlopI feOOlbf ). 


PG1307+085, z = 0.155 

Th e host galaxy of this radio quiet quasa r is an elliptical wit h no evidence of a disturbance (B ahcall et all 
Il997t iHamilton et al.1 [2002: Gu von et, all 120061 IVeilleux et al.l I2009I ). The center of the z = 0.14 cluster GMBCG 
J197.43444+08.33445 dKoester et al.l I2007D lies at a projected separation of only 43.5" from the quasar. Given the 
redshift separation it is unlikely that PG1307+085 is a cluster member but many galaxies in the 90Prime field likely 
are. 

Our results for B gq in the r band are consistent with IMcLure fc Dunlopl ( 2001b ) and I Finn et al.l ( 2001 ) within the 
large uncertainties. The i band result for B gq for this field is very large, 1461+778. The B gq value found for the g 
band, 259+120, is also large relative to other filters. The presence of GMBCG J197.43444+08.33445 in the foreground 
may be influencing the galaxy counts for this field, though it is not clear why this should only be the case for the g 
and i filters and not u and r, which both show B gq consistent with B gg . 


PG1404+226, z = 0.098, NLS1 

Hoi (I2005I) conclude that the host g alaxy of this quasar has a very low star formation rate, <0.14 M 0 yr 1 , based 
on the [O II] measurement of lKuras zkiewicz et al.l (I 2000D. It also h as a low black hole mass. log(MBH/Mg) = 6.88, as 
estimated from the H/3 line width ( Vestergaard fe Petersonll2006l) . iMcLeod fc Riekel (Il994al) find that an exponential 
profile with a scale length of 2.98 ft” 1 kpc gives a reasonable fit to their H band image of this AGN. 

There a re ~5 0 clusters, gro u ps, or candidates identified in the 50' around PG1404+226, including five compact 
groups of iMcConnachie et ahl (20091) with separations less than ~22' and unknown redshift. The cluster WHL 
J140609.9+221605 has a separation of 8.17' from the quasar and a redshift of z = 0.393. 

All values of B gq for this field are consistent with B gg , except in the g band, where we find B sg =942±522. 


PG1444+407, z = 0.267 

The morphology of the host galaxy of this radio q u iet AGN is so mewhat ambiguous. iBahcall et al.l (19971) confirm 
the possible presence of a bar sugge sted bvlHutchings fc Ne ff (1992) as the result of an old merger stage. Interestingly, 
based on the [O 11] measurement of lBaldwin et al.1 (~1989| ). iHol (120051 ) infer a fairly high star formation rate in the host 
of 19.4 Mg yr ”, on the order of that found for interacting and luminous IR galaxies. 

There are 48 galaxy clusters, groups, or candidates within 50' of the quasar, with the closest, ZwCl 1443.8+4043, 
lying at a separation of 12.6', but unknown redshift. 

For this field we find B gq values in the g , r, and i filters consistent with B gg . T he r band resu lt i s consi stent 
with small posi tive values with large statistical uncertainties from the other studies of lMcLure fc Dunlopl ( 2001bl ) and 
iFinn et ahl (20011) . In the u band, we find a B gq B gg by more than 2cx, reflecting the general trend for the u band 
to show the smallest values of B gq , though we note that the limiting u magnitude for this field corresponds to M*+1.4 
at the redshift of this AGN rather than to M*+2. 


PG1545+210 (3C 323.1), z = 0.264 

This FRII radio loud AGN has an elliptical host galaxy (McLeod fc Riekel 1 1994bl : IBahcall et al.l[l997l) with a faint 
elliptical companion galaxy at a projected separation of 7.7 ft” 1 kpc. IHoI ( 2005I) finds a moderately large star formation 
rate of 9.7 M 0 yr” 1 in the host galaxy, based on the [O III measurement o f IWi lls et all ( 1993 ). From population 
synthesis model fits to the optical spectrum of the companion. iCanalizo fc Stocktonl (I1997D report, evidenc e of a possible 
intera ction 2.3 Gyr ago. This qua sar was found to have a rich cluster environment (Abell class 1) by lOemler et al.1 
(fl972l) . and iKirhakos et al.l (1994D find 22 likely galaxies in their optical images described above in Sec. IA.1I It is 
located 1.96' away from another broad line AGN at a very similar redshift (SDSS J154749.70+205056. 5, z = 0.26547). 

These AGNs are likely members of the X-ray indentified clu ster ZwCl 1545.1+2104 (z = 0.266, iPiffaretti et al.1 
(2011D ). also found in the SDSS DR7 data (z = 0.26735. lHao et ahl (2010D ). There are ~70 other galaxy clusters, 
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groups, or candidates within 50' of PG1545+210. We note that this field also contains the highest number of 2 < 0.5 
AGN of any of the other 90Prime fields, listed in Tabic [9] 

Here we a B gg value in g consistent with B gg , and the r band value identically zero, as for MRK 586. This is in 
agreement wit h the HS T/WF PC2 F606W results of lMcLure fc Dunlod : 2001 b :. However, given the smaller uncertainty 
quoted by i Finn et al.l 1|2001| ). our value r band value differs from their positive results by ~ 3cr. As for PG1307+085, 
we find a very large value of B gg for this held in the i band and a much smaller but significant signal in u, consistent 
with its known cluster environment. The limiting u magnitude for this held corresponds to M*+1.7 at the redshift of 
this AGN, rather than reaching to M*+2. 


PG1612+261, z = 0.131 

iMcLeod fc R ickc (1994a) hnd that an exponential prohle with a scale length of 1.16 h~ 1 kpc gives a reasonabl e ht 
to their H band image of this AGN and its im mediate environment meets the criteria of lMcLeod fc Riekel (11994H ) for 
bright interacting companion. (See Sec. IA.3I 1 

It may be associated with the cluster NSC J161419+260832 {z = 0.1396) but given the velocity separation and the 
^0.5 /i -1 Mpc separation on the sky from the reported cluster center, it is likely on the periphery. There are ~47 
other unique galaxy clusters, groups or candidates within 50' of the quasar found in NED. 

The g band B gg value for this field, 354±237, is marginally greater than B gg , but it is consistent with B gg in u, r 
and i. The most closely separated z < 1 QSO in this field, SDSS J161335.33+263127.7, lies at z = 0.708 and ~28.5'. 

[HB89] 2U1+175, z = 0.211 

Th i s radio loud AGN has an elliptical host, with some suggestion of an interacting companion dHutchings fe Neffl 
Il992t iDunlon et al.ll2003l h There is only one known galaxy cluster in this field, within 46' of this quasar. _ 

We h ave no i band image from 90Prime to compare with the flS'T/WFPC2 F675W result of Mc Lure fc Dunloul 
(I2001bf ). but our large r band result, 991+355 is larger than theirs by >2.5 a. 
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